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An Apparatus for Measuring Thermal Expansion at 
Elevated Temperatures 


B. D. Rothrock and R. K. Kirby 


Institute for Basic Standards, National Bureau of Standards, Washington, D.C. 20234 


(November 23, 1966) 


An apparatus has been developed for making precise measurements of the absolute thermal 


expansion of refractory materials at temperatures up to 1600 °C. 
comparator and a controlled gradient vacuum furnace. 


This apparatus consists of an optical 
Special care was taken on the techniques 


of measuring length and temperature and on the calibration of the microscopes and thermometers. 
The expansion data obtained with this apparatus on specimens of platinum and sapphire have a stand- 
ard deviation of less than 30 ppm and are estimated to be accurate to within 50 ppm. 


Key Words: Controlled-gradient vacuum furnace, elevated temperatures, optical comparator, 


thermal expansion. 


1. Introduction 


Accurate data on the thermal expansion of solids at 
temperatures above 1000 °C are needed for several 
reasons. (1) They are needed on at least one repro- 
ducible material so that engineering data on other 
materials can be obtained with relatively simple 
apparatus by a comparison method [1].! (2) The 
material on which the data is obtained can be used to 
measure temperatures in high-temperature x-ray dif- 
fraction cameras and diffractometers [2,3]. (3) Accu- 
rate high-temperature data can be useful in the study 
of lattice defect formation [4], phase transitions [5], 
and the theory of thermal expansion [4, 6]. The appa- 
ratus described in this paper was constructed primarily 
to furnish accurate thermal expansion data up to 
1600 °C on a reference material that could be used by 
others in measuring thermal expansion. 
this, special care was taken on the techniques of meas- 
uring length and temperature and on the construction 
of the furnace. 

An optical comparator method was chosen because 
it can be used to make accurate and precise measure- 
ments of the expansion of solid material at elevated 
temperatures. The present method makes use of a 
pair of microscopes mounted on an invar bar and each 
measurement of the length of the specimen is com- 
pared to the length of a reference bar kept at room 
temperature. This procedure gives a correction that 
compensates for any apparent change in distance 


between the microscopes caused by the expansion of 


' Figures in brackets indicate the literature references at the end of this paper. 


Because of 


the invar bar or the tilting of the microscopes. The 
nominal length of the specimen between the fiducial 
marks is 10 cm. The specimen, hung vertically in 
the furnace, is free of all restraint so that measure- 
ments of expansion will be accurate as long as the 
specimen does not creep significantly under its own 
weight. The vacuum furnace is specially designed 
to correct for end heat losses. In addition the tem- 
perature profile over the length of the specimen was 
measured with a calibrated probe thermocouple. -A 
detailed description of the apparatus and techniques 
is given in the following sections. 


2. Apparatus 


2.1. Furnace 


To make accurate length versus temperature meas- 
urements the furnace must be designed and con- 
structed to produce a uniform temperature along the 
length of the specimen or the temperature gradient 
along the specimen length must be known at each tem- 
perature to determine the effective or average tempera- 
ture of the specimen [7]. To obtain a_ uniform 
temperature over the specimen length a controlled- 
gradient furnace was constructed that has five indi- 
vidually controlled heating sections obtained by 
tapping an evenly wound resistance heater at 5 cm 
intervals. By the use of isolation transformers the 
power in any one section can be controlled with only a 
negligible change in power in the remaining sections. 
This design, as shown in figure 1, makes it possible 
to adjust the power input at any temperature to com- 
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FiGuRE 1. Schematic diagram of heater control circuit. 


pensate for end heat losses. A voltage regulator was 
used to compensate for line voltage variations as it 
was felt that this would be the main cause of tempera- 
ture variations. The power required to obtain 1600 °C 
was 2.6 kW. 

The heater wire, having a diameter of about 1 mm, 
is wound about a core in a 1.6 mm deep continuous 
groove, 6 revolutions per 25.4 mm. Tantalum, tungs- 
ten, and rhenium wire have been used. Although 
each have desirable resistivity and high temperature 
properties, rhenium proved to be the best because of 
its high-temperature inertness and greater ability 
to be worked. 

The high-purity alumina (99% Al»O3) core on which 
the heater is wound is 30.5 cm long and has a 1.91 
em LD. and a 2.54 em O.D. Sight holes, 6.3 mm in 
diameter, in the core, see figure 2, are alined with 
holes in the heat shields and the windows in the outer 
shell. Two windows on one side, with centers 10 
cm apart, allow viewing of the specimen for length 
measurements. A third window on the opposite side 
permits temperature measurements with an optical 
pyrometer. 

The specimen, a 12 to 14-cm rod having a diameter 
of about 5 mm, is hung from the top of the furnace by 
two loops of tungsten wire. The small lower loop 
runs through a small hole in the top of the specimen 
and is connected to the larger upper loop by a tungsten 
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joint which allows the specimen to hang freely in a 
vertical position. Vertical, lateral, and_ rotational 
movement can be achieved through a bellows and 
vacuum seal. 

The control and measuring thermocouples are po- 
sitioned in 0.76 mm holes in the wall of a high-purity 
alumina tube. This tube, 1.11 em [.D. by 1.75 cm 
O.D., is 35.6 cm in length. The thermocouples are 
shielded by a grounded rhenium foil placed between 
the heating core and thermocouple protection tube. 
A small a-c pickup can be detected but did not affect 
the thermocouple emf output as read by either the 
potentiometer or recorder. 

The outer stainless steel shell of the furnace is 
water-cooled and vacuum sealed. The vacuum system 
consists of a fore pump, oil diffusion pump, water- 
cooled baffle, and vacuum gages. The vacuum gage 
indicator, power supply and controls, and temperature 
indicators are mounted in a relay rack. 


2.2. Comparator and Reference Scale 


The comparator shown in figure 3, consisting of 
two microscopes mounted on an invar bar, was ob- 
tained commercially. The microscopes have filar 
micrometer eyepieces, 50X magnification, and relay 
lenses which allow a working distance of l6cm. They 
can be raised or lowered and rotated about a vertical 
axis. The comparator, however, was modified to 
allow for precise alinement after being rotated between 
the specimen and scale positions. 

The reference scale, mounted outside the furnace, 
is made by wrapping a 0.020 mm tungsten wire about 
a vitreous-silica rectangular bar 2.54 cm xX 1.91 
cm X15.2 em. The wires are spaced 1 mm apart 
by grooves at the edges of the silica bar and held in 
position by epoxy. The back surface of the clear bar 
was ground so that when illuminated from the back 
the wires stood out sharply against an evenly illu- 
minated field. The scale was calibrated accurately to 
within +0.2 micrometers by using a line standard, 
Decimeter No. 43, provided by the Length Section 
of the National Bureau of Standards. The position 
of the scale can be adjusted in the vertical and hori- 
zontal directions. 

The scale eliminates any apparent or real shift in 
the measurement of length caused by expansion of 
the invar bar or tilting of the microscopes and estab- 
lishes the absolute length of the specimen. In pre- 
liminary tests on the microscopes, measurements were 
made on the reference scale while changing the 
ambient room temperature. These tests disclosed 
that one of the microscopes was tilting in its mount 
and might have contributed an error in length measure- 
ments of about 100 micrometers/meter. The mount 
was remachined and provided no further difficulty. 
This experience emphasizes the usefulness of a refer- 
ence scale for if the defective mount had not been 
detected the use of the scale would have automatically 
corrected the length measurements. 

The microscopes were alined and calibrated by 
first adjusting the invar bar on which the microscopes 
are mounted to the vertical position. Next, the micro- 
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K. Alumina core 


L. Specimen 
M. Window 

N. Shutter 

O. Vacuum port 


P. Vicor tube 

Q. Vacuum 

R. Thermocouple protection tube 
S. Pyrometer viewing hole 

T. Tantalum tube. 





FIGURE 3. Photograph showing relative position of comparator, 


furnace and reference scale. 


scope tubes were adjusted until they were horizontal 
and parallel. The microscopes were then moved in 
their individual mounts until they were both focused 
on a freely hanging 0.38 mm wire. The paired micro- 
scopes were turned toward the reference scale and 
the position of the scale was adjusted until it was 
focused in each microscope. With the microscopes 
focused on the scale each eyepiece was calibrated by 
measuring the separation of two wires. After many 
repeated measurements the values for the eyepieces 
were determined to be 1.128 and 1.095 micrometers 
per division. 


3. Control and Measurement Techniques 


3.1. Temperature 


To control the temperature gradient in the furnace 
four differential W-Re thermocouples were placed in 
the thermocouple protection tube. One of the junc- 
tions of each of the thermocouples was placed at the 
center of the furnace and the other junction was 
placed at the center of one of the other four heating 
sections. This provides a comparison of temperature 
at the center of each heating section with the tempera- 
ture at the center of the furnace. 

A multipoint strip chart recorder is used to record 
the output of the differential thermocouples. The 


recorder scale is zero centered with a full scale deflec- 
tion of +300 uwV (10 wV per scale division). The plot 
can be read to better than 1 scale division and hence 
the differential temperature to better than 1 °C. The 
24-second printing cycle of the recorder is about the 
response time of the furnace to a power change. 
Using this system it is not difficult to maintain a 
gradient within +2 °C as indicated by the recorder. 
In the time it takes to make a set of readings the 
average temperature usually did not vary more than 
=3°C. 

To test the actual temperature profile at the position 
occupied by the specimen an arrangement was made 
for measuring the temperature with a single absolute 
thermocouple. This probe thermocouple was inserted 
into the furnace through the top plate by means of a 
sliding vacuum seal, see figure 2. The wires were 
insulated with a two-hole high-purity alumina protec- 
tion tube. The junction was about 5 mm below the 
end of the alumina tube but was shielded by a thin- 
walled 4-mm O.D. tantalum tube. The tantalum tube 
extended about 3 c.1 above and below the junction 
and had two small holes in it to allow the optical 
pyrometer to “see” the junction when it was located 
at the center of the furnace and to use the tube as a 
blackbody when the junction was located slightly above 
the center of the furnace. 

The probe thermocouple measured the temperature 
profile along the axis of the furnace at several different 
temperatures. At each temperature the current in 
each heating section was adjusted until the print out 
on the strip chart recorder indicated that all four dif- 
ferential thermocouples were at the same temperature. 
The probe thermocouple was then raised and lowered 
and readings taken at 1 cm intervals. As shown in 
figure 4, the actual temperature profile has an overall 
slope with dips at the viewing windows. This same 
slope was found when the probe thermocouple was 
inserted into the furnace through the bottom plate. 
The temperature profile that is obtained when each 
heating section has the same current flowing through 
its heating element is also shown in figure 4. 

Temperature measurements are made at the center 
of the furnace with a W-Re thermocouple and a Pt-Rh 
30-6 thermocouple placed in the thermocouple pro- 
tection tube. These thermocouples and the optical 
pyrometer were calibrated by comparison with a probe 
thermocouple that had been calibrated by the Tem- 
perature Physics Section of the National Bureau of 
Standards. During this calibration the probe thermo- 
couple was at a point 2 cm above the center of the 
furnace where the slope of the temperature profile 
was nearly zero. The value of the average tempera- 
ture and its position was determined from the tem- 
perature profile over the specimen length. A correction 
was made to the indication of the probe thermocouple 
corresponding to the difference between the average 
temperature and the temperature at the 2 cm position 
as taken from the temperature profile. 

Because the thermocouples were calibrated with 
the probe thermocouple occupying the same space 
that the specimen occupies during an expansion test, 
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FIGURE 4. Temperature profile along the axis of the furnace (top 
of furnace to right): A, current adjusted in each heating section 
until all four differential thermocouples were at the same tempera- 
ture; B, same current flowing through each heating section. 
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represent 


the temperature 
the 


it is assumed 
accurately 
specimen. 

While the temperature measurements made with 
the optical pyrometer and the thermocouples were 
usually in reasonable agreement, within +3 °C, the 
measurements made with the pyrometer were only 
used as an additional check on the thermoucouples. 
After the furnace had been heated above 1500 °C, 
slight deposits could be seen on the windows and the 
temperature measurements made with the pyrometer 
would be as much as 50 °C lower than those made 
with the’ thermocouples. 


measurements 
temperature within the 


3.2. Length 


In order to obtain reliable measurements of length 
the fiducial marks should be an integral part of the 
specimen or not subject to arbitrary shifts. In order 
to obtain measurements of high precision the fiducial 
marks should be small (of the order of 0.025 mm when 
using a magnification of 50X), symmetrical, and 
capable of forming an image with sharp focus. Sev- 
eral experimenters [8, 9, 10] have used fiducial marks 
made by peening small wires into holes drilled into 
the specimen, clipping the wires off close to the 
specimen and placing nicks in the resulting sharp 
edge. Other investigators have used marks scribed 
or indented into the surface of the specimen [4] or 
sighted on the ends of the specimen [11, 12]. 


89 





























(c) 








ie lad 


FIGURE 5. Fiducial marks: (a) 0.020 mm tungsten wire tied taut 


around the specimen; (b) platinum specimen; (c) sapphire specimen. 


A promising method of tying two 0.020 mm tungsten 


wires taut around the specimen was tried, see figure 
5 (a). Each wire crossed in front of a 1.5 mm diam 
hole drilled in the specimen and was held in position 
by running through grooves on either side of the hole. 
When the furnace became hot enough that the wires 
were self illuminating they stood out sharply against 
the dark background and precise settings, within 
+0.7 micrometers, with the filar eyepiece could be 
made. At temperatures above 1500 °C, however, 
probably because of recrystallization, the wires grew 
in length and warped. Because they were no longer 
straight any further readings could not be relied upon. 

It has been found that small holes (of the order of 
0.25 mm) drilled through a specimen [13] or small 
nicks or grooves [14] will serve as reliable fiducial 
marks. It is sometimes difficult, however, to obtain 
a round hole that is small enough for high precision 
measurements. 

Satisfactory fiducial marks were finally made by 
placing a small symmetrical nick at the edge of the 
specimen. This was done for the platinum specimen, 
see figure 5 (b), by drilling a 1.73 mm hole radially 
into the rod to within about 0.25 mm of penetrating 
a milled flat on the opposite surface. A second hole 
of 1.57 mm diam, concentric with the first hole, was 
drilled in the flat, penetrating through to the larger 
diameter hole. This left a thin ring around the inside 
of the hole and a small nick was placed in the ring 
with a hardened steel ruling edge to serve as a fiducial 
mark. In the case of the sapphire specimen, see 





figure 5 (c), a local area on one side was ground with 
a diamond wheel until a thin ridge was left, the faces 
of which were parallel to the axis of the specimen. A 
small nick was then placed in the ridge with a 0.25 
mm carbide disk as a fiducial mark. The precision 
of setting the cross hairs on these marks was found 
to be comparable to that of setting on the tungsten 
wires. 


4. Procedure 


With the specimen in the furnace and the furnace 
sealed the vacuum pumps were turned on and the 
pressure reduced to about 1 X 10-5 torr. Power was 
applied slowly to heat the specimen. When the de- 
sired temperature was reached the optimum tempera- 
ture gradient was obtained by adjusting the power sup- 
plied to the two heating sections at either end. This 
condition was held for 10 min to allow the specimen to 
stabilize. A longer stabilizing period was not neces- 
sary since it was found that no change in length of the 
specimen could be detected after a 5 min wait. The 
microscopes were focused on the specimen by adjust- 
ing the position of the comparator and on the reference 
scale by adjusting its position. The emf of the two 
absolute thermocouples was determined with a poten- 
tiometer to the nearest microvolt and the brightness 
temperature of the specimen was measured using the 
disappearing filament optical pyrometer. Five set- 
tings with each microscope were made in turn on the 
specimen, the reference scale, and again on the speci- 
men. Finally, the temperature measurements were 
repeated. 


5. Discussion 


Although care has been exercised in measurement 
of temperature and length, errors inherent in any 
system operating at these temperatures restrict the 
accuracy in the measurement of these values and the 
determination of the coefficient of expansion. An error 
in temperature measurement of 1 °C in a specimen 
having a coefficient of 10 X 10-§/°C is equivalent to 
an error in length measurement of 10 micrometers/ 
meter. Also, an error in temperature measurement 
of 1 °C, or in length measurement of one micrometer 
in a specimen 10 cm long having a coefficient of 
10 x 10-§/°C and heated over a temperature range of 
100 °C, will result in an error of less than +0.1 X 10-8/ 
°C in the determination of the coefficient. 

Uncertainties in temperature measurements may 
arise from a shift in calibration of the thermocouples 
through contamination and aging. Although high- 
purity alumina (99% Al,O 3) has been used fer thermo- 
couple insulation the trace impurities can affectively 
change the output of the thermocouples. To reduce 
the effect of this and other possible sources of error 
the following procedure was used: (1) The probe 
thermocouple was calibrated by the Temperature 
Physics Section of the National Bureau of Standards, 
(2) the absolute thermocouples were calibrated by 
comparison with the probe thermocouple, (3) the 


expansion tests were conducted, (4) the absolute 
thermocouples were again calibrated, and (5) the probe 
thermocouple was again calibrated. This procedure 
demonstrated that the uncertainties in temperature 
measured during the thermal expansion tests did not 
exceed 2 deg at 800 °C and 4 deg at 1600 °C. The 
standard deviation of the differences between the 
temperatures indicated by the two absolute thermo- 
couples was 1.5 °C, and that of the average of the two 
readings is 0.7 °C. 

Preliminary results? have been obtained on speci- 
mens of platinum and sapphire. Figure 6 indicates 
the difference (residuals) between observed expan- 
sion of these materials and the expansion calculated 
from quadratic equations fitted to the data by the 
method of least squares. For the platinum and sap- 
phire specimens the standard deviations computed 
from the residuals are 23 and 27 micrometers/meter 
respectively. 

From these estimates of the precision of the tem- 
perature measurements, + 0.7 °C, and of the expansion 
measurements, +23 and +27 micrometers/meter, the 
precision of the length measurements expressed as 
standard deviations are 2.1 and 2.6 micrometers, 
respectively. . 
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FIGURE 6. A plot of the differences between observed expansion and 
the expansion calculated from an equation fitted to the data by the 
method of least squares. 


Assuming that the accuracy and precision of the 
length measurements are comparable and the uncer- 
tainties in the temperature measurements are as 
previously stated, thermal expansion measurements 
made with this equipment can be estimated to be 


2 A complete report on these measurements will be made after the investigation has been 
completed. 
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accurate within +30 micrometers/meter at 800 °C 
and within +50 micrometers/meter at 1600 °C. 

It was discovered that a small wedge (6 min) be- 
tween the faces of each vitreous silica window con- 
tributed an error in the measurement of expansion 
from room temperature. Since the measurement of 
initial length is made with the specimen outside the 
furnace there is no precise measurement of the change 
in length of the specimen from room temperature to 
900 °C where the specimen is self illuminating. This 
error was virtually constant for a fixed position of the 
window, however, and did not significantly affect the 
determination of the coefficient of expansion at the 
operating temperature. To eliminate this and other 
possible errors two higher quality windows are being 
made that will have surfaces flat to 1 wavelength of 
visible light and parallel to 1 sec. 

Light sources are also being placed in the furnace 
so that the length of the specimen can be determined 
inside the furnace at ambient and intermediate 
temperatures. 

With the improvements discussed in the previous 
paragraphs and the preliminary results obtained on 
platinum and sapphire it is believed that the apparatus 
and the procedure described in this paper will yield 
important and useful thermal expansion data on high 
temperature materials. 


The authors thank L. Frenkel for contributions made 
in the design of the apparatus and to T. A. Hahn for 
helpful discussions. 
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Determination and Smoothing of Fourier Coefficients 
Representing Piecewise Continuous Functions 


B. A. Peavy 


Institute for Applied Technology, National Bureau of Standards, Washington, D.C. 20234 
(February 7, 1967) 


This paper presents a method of solving for Fourier coefficients where the dependent variable 
can be expressed as a piecewise continuous function, when various conditions of continuity and 
smoothing are assumed. An example is included to show the effect of smoothing in the region of a 
discontinuity for a system composed of two materials that exhibit a discontinuity at their interface 
and surrounded by a third material which does not have a discontinuity. An advantage to be gained 
from smoothing is an increase in the convergence of a finite Fourier series representation of a piece- 
wise continuous function in the region of the discontinuities. 


Key Words: Continuous function, Fourier series, piecewise, smoothing. 


1. Introduction 
Many time-independent problems of a practical nature involve the solution to the partial 
differential equation 


02v . O7v 


“i, 
+—=0 


ox? AF Oz 
in rectangular coordinates, or 


O*v l dv l Ov Pv 
eg ae ene 


or ror roe a2 
in cylindrical coordinates. From the general solution, one particular solution is in the form 


: . nz 
v=B+Cz+$) An sin : (2) 
n 


@ 


where the Fourier coefficients, 4, =A,(x, y) for rectangular coordinates, or A, =A,(r, ¢) for cylin- 
drical coordinates, are determined from boundary conditions in the region 0 <z<@, and fixed 
coordinate positions for x and y, or r and ¢g, with the condition that other necessary solutions of 
(1) disappear on this boundary. Also, B+ Cz must satisfy v at the end points, z=0 and z=wo. 

From an experimental point of view, the dependent variable, v, can be approximated from 
measured values at several z positions on the boundary by a relationship, v=g(z). The Fourier 
coefhicients are usually simply determined if the function g{z) and its derivatives are continuous 
over the region 0 <z<q@. For various applications, the variable v must be expressed as piecewise 
continuous functions in this region, with apparent discontinuities at discrete points. If it is known 
that the discontinuities do not exist, then an appropriate smoothing function can be employed. 

The purpose of this paper is to solve for Fourier coefficients, when it is assumed that the 
dependent variable v can be expressed by polynomials representing discrete portions of the region 
0 <z<a@, and that conditions of continuity and smoothing exist. An application of functions 
developed to a hypothetical problem is also presented. 
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2. Analysis 


On a fixed surface lying in the region 0 < z 


v=gi(z): 


a, 3, 


have discontinuities at the points z;, i= 1, 2, 


where i=l, 


. m, z=0, and zm=w. 


<w, the dependent variable v is defined 


Zi—1 25 Zs 
The functions gi(z) and their derivatives may 


m—1l1. Substituting in (2), multiplying through 


by sin kzrz/w, and integrating between the specified limits gives the following equations 


i-1 


Adding the above equations 


> naz . #82 
SiAn sin sil z=— 
n 0 @ @ 0 


from which 


9 9 
An=— {Cw(—1)"—B[1—(—1)"]}+—- 
nt ws 


Repeated integration by parts yields 


” (B+Cz) sin aan 


kaz 


2 _ naz, 
[" te ‘gi(z) —B—Cz] sin“ qe= ya. |" sin ——~ sin dz. 


“i-1 


m 25 koz 
= TZ 
gi(z) sin — dz, 
W 


ad 


(63) 2j_1 


| nz 
gi(z) sin dz. 
W@W 





. NZ 
[ee sin 2 


pve 


where 
Letting 


superscripts ' 


NZ; 
sin 
@ 
Gi =—— 4.2. (zi) — 2, (zi — 
("2)’ 6i Bi+1\4i 
@ 


@ 
a {gz —giee)— 


g; (0) 


, represent the first, second, third, et 





‘., derivative with respect to z. 


[gi (zi) — 23 (zi) 





- ge 


ie). 
@ 
[gi'(zi) — git (zi) 


a] 





gi'(0) 


Ge) 

@ @ 

_ Bm(@) gino) 

(= F =) 
@ @ 
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€ 


2 m—-1 
{g(0)— B+(—1)"[Cw+ B— gn(w)|— H+(-— 1)"K}+ an pa Gi. 


i=1 


A, 


nt 
For the condition that B= g,(0) and Cw + B=g,(@), the Fourier coefficients are 


> 2 m-l ' 
A,=— {(—1"K—H}+-— eA Gi. 
nit @ 


i=1 


Assume for the region, 0 < z S , that v as a function of z behaves as shown in figure 1, where 
£1, 82, 83, and g; are known functions of v, (perhaps as determined from experiment) and fi, f2, and 
fs encompass regions of uncertainty or regions which cannot be defined by experiment. If it is 
known that there are no discontinuities within the regions of the uncertainties, it is possible to 


smooth or bridge between the known functions with a cubic polynomial which gives continuity to 


two adjacent g functions and their first derivatives. Examining the first region that includes an 
uncertainty, the cubic polynomial, f:(z)= Mo + Miz+ M2z? + M3;z' is solved for the coefficients, M, 
such that gi(z1)=filzi), gi(zi) =Sf1 (21). go(z2)=fi(ze) and gi(z2)=/{ (zz), from which the necessary co- 
efficients are 





“—— 3[ge(z2) — gi(z1) Jz — €[(3Z + ©)gi (z1) + (37 — Ege (z2)] 
a de 





€[gi(21) + 8(z2)| +21(21) —go(z2) 
4e3 


M, ———— 


where 27=z; +€=22—€ is the midpoint between z; and z.. Assuming bridging between g; and ge, 
then 


G — F, + Fy 


where 





FIGURE 1. Variation of v with respect to z showing the 
defined functions g). go. g3, and g4 connected by functions 
f,. fo, and fs through regions of uncertainty. 

















Making the substitutions z;=z—e€, z.=z+e, f/'(z) =2M2+6Msz, and f/''(z) =6Ms in F»2, yields 


_(nie\ . 
Ul (=) sin 


a , ’ ‘ Ww 
F.= [gi (z:) —gs(z2)] —————_ 


ei] 





— [g1(z1) —ge(z2) + e{g}(z1) + 23(z2) }] 
where 


u(p) =“BP=1-F 4 


T(B) = 33 (sin B—B cos B) 





Plots of U(B) and T(f) are given in figure 2. 

For the interval z; < zz, (4b) is the bridging or smoothing function that gives continuity 
to gi(z) and g2(z) and their first derivatives at z; and z2, respectively. This function also eliminates 
computation of the cubic polynomial f;(z), where the values of M2 and M3; sometimes give inordi- 
nately large values with subsequent loss of significance. The same analysis can be applied to the 
regions encompassed by f2(z) and f3(z) so that G2 and G; can be found. 


3. Specific Cases 
Where successive differentiation terminates, the functions giz) may be represented by poly- 
nomials. There are many variations of gz) possible, but probably the best example is to examine 
a set of straight lines as given in the following specific cases. 


3.1. Functions as Straight Lines 


The functions gi(z) are represented by a set of straight lines. Here g/'(z)=H=K=0, so that 
(3) becomes 


nS, 
An=— > Gi. 
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FIGURE 2. Variation of U(B) and T(B) with B. 


FIGURE 2. Variation of U(B) and T(B) with B. 


If a set of straight lines which do not intersect at critical positions z;, as shown in figure 3 are 
assumed, then 


W 


Gi= (—] [27 (zi) — gj+(zi)] sin 
nt 


N77Zj @ NWTZ;j 
= tele my, dadloee. ee (6) 
@ nit 63) 


where the functions and derivatives are evaluated at z;. 


For the cases where the lines do intersect 
ee ee = pi 


. i= gi+i, the second term on the right side of (6) is zero. 


3.2. Straight Line Functions Smoothed Near Intersections 


The convergence of (6) when substituted in (5) and (2) is not very satisfactory, and it is possible 


to increase the convergence by joining the straight lines with smoothing curves that give continuity 


to the lines and their first derivatives at points z;— y; and z;+€;, as shown in figure 4. 
Proceeding as in the derivation of (4), where F', =0, yields 


T(nXj) COS Qj 


9 ("=) 
@ / 


Gi = [(€i — yi) (gi. (zi) — Bi(zi)) + Agis (zi) — gilzi))] 


U(nXj) sin aj 


9 
(=) 
Ww 


— {gi,1(zi) — gi(zi)} 


where 


As €;=yi > 0, (7) becomes (6). 


There are many possible specific cases for (7), i.e., €¢;=0, 
yi=0, €:=yi, ete. 


Essentially, smoothing increases convergence by introducing the factors T(B) 
and U(8), 8=ndAj, to the terms on the right side of (6). The plots 7(8) and U() (fig. 2) indigate that 
the advantage of the increase in convergence is not realized until B > 0. 
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FIGURE 3. Variation of v as a set of straight lines with FIGURE 4, 


Variation of v as a set of straight lines with 
function discontinuities at z, and 22. 


smoothing in the regions of 2, and 2». 


3.3. Example of Smoothing 


Figure 5 is used as an example of the variation of the dependent variable v with the dimension 
z in order to show the need for smoothing in the region of a discontinuity. For a hypothetical 
process, consider the potential on the surface of a system of two materials which is surrounded 
by a third material (not a perfect insulator), and where, on the surface of the system, the region 
0 <z=<0.5 represents the potential in a perfect conductor, the region 0.5 < z S 1.0 represents that 
in a material of finite conductivity, and the discontinuity represents the contact resistance between 
the two materials. Usually, experimental processes allow measurement of potential in the two 
regions, but measurements at or near the discontinuity are not always possible or desirable. As- 
suming homogeneous isotropic materials, measurements taken in the two separate regions give 
linear relations which, when extrapolated to the point of contact (z=0.5), do not meet in a common 
point and give an apparent discontinuity. When it is known that the process must be continuous, 
i.e., the surrounding material is continuous at z=0.5, then a smoothing function may be applied. 

Numerical values for substitution in (7) are z;=0.5, g:=1.0, gj =0, ge =0.95, g3=1.9, and 
€=y, where a=n7/2, \=7re. Substitution in (7), (5) and (2) gives 





y sin 27nz _ 3.8U|(2n— sin (2n — l)zz 
y= 1-2-5 (— 1p | sn 27n _ 3.801 n— 1)A] sin (2n—1)77z). (3) 


20n7 (2n — 1)?77r? | 


Figure 6 is a plot of v in the region of the discontinuity shown in figure 5 for the summation of 200 
terms of the above series and values of €= 0.0, 0.002, 0.003, and 0.004, and their respective smooth- 
ing curves. Figure 6 shows the influence of the increased convergence of the series toward its 





1.0 





FIGURE 5. Variation of v with z for the problem of 
section 3.3. 











FIGURE 6. Numerical values of v from (8) in the region 
of the discontinuity of figure 5 for a summation of 200 
terms and values of €=0, 0.002, 0.003 and 0.004. 
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intended values with increase in the value of €. For this case, fairly good convergence of the series 
is attained for n > 1/2e. For the summation of 300 terms, the fluctuations shown on figure 6 are 
virtually eliminated, except for «= 0. 


4. Discussion 
The example of section 3.3 as shown in figure 5 was chosen because it illustrates an extreme 
condition that a physical process may undergo: that is, the dependent variable has a zero slope 
up to the discontinuity at the point z=0.5. It is then to be expected that actual physical processes 
where apparent discontinuities occur would give better convergence when smoothing conditions 
are employed. 


Generally, the Fourier sine series (2) can be used to represent a piecewise continuous function 


over an interval, 0 <z<@, by the use of (3). Numerical evaluation will, of course, present prob- 
lems of convergence in the regions where discontinuities occur (as shown in fig. 6, €=0). The 
lack of sufficient convergence in (2) can cause serious problems when (2) is used in the solution 
of (1). It is therefore expedient to employ the bridging or smoothing relationship (4) (or (7) for 
the special case of a set of straight lines) for which adequate convergence can be assured for a 
reasonable number of terms in the series of (2). 


A fairly simple heat transfer system that may be analyzed employing the method presented 
in this paper is shown in figure 7a by a cross section through a composite cylinder. The inner 
cylinder (regions 1 and 2) 0 <r <a is composed of two materials of different thermal conductivities, 
usually with values greater than that for the material of the outer cylinder (region 3). For the 
determination of thermal conductivity, region 1 may comprise a material of known thermal con- 
ductivity and region 2 a material of thermal conductivity to be determined. It is necessary to 
determine the heat transfer characteristics in the material of region 3 in order to evaluate the radial 
heat loss or gain of the inner cylinder, 0<rSa. 

The temperature distribution at r=a will vary according to the thermal conductivities of the 
materials of regions 1 and 2, and the distribution may have an infinite number of variations. One 
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FIGURE 7. Heat transfer system. 


a) Cross-section through composite cylinder 


b) A poss ible variation in temperature potentia alr aandr 


such distribution for r=a with a thermal contact resistance at z=0.5 is shown in figure 7b, where 


the intended heat flow is in the positive z-direction with a unit temperature potential at z=0 and a 


zero value atz=1. On r=b the temperature is assumed to decrease linearly in the above limits. 
9 


considering region 3 only. a solution of (1)(assuming no @ dependence) is 


ee. AF (nar, nzb) sin n7z 


— 





Fy(nza, nb) (9) 
where Fo(x, y)=I1o(x)Koly) — Io(y)Ko(x) and I) and Ky are modified Bessel functions of zero order and 
first and second kind, respectively. The substitution of r=a in (9) gives a form similar to (2), for 
which a solution for 4, can be found from (7), the temperature distribution shown in figure 7b and 
a reasonable value of € that insures convergence of the resulting series. 

Two of the objectives of the NBS Heat Transfer Laboratory are to develop and analyze methods 
for the determination of thermal conductivity of materials which vary widely in thermal properties 
and to make measurements of thermal conductivity for government agencies. For both analysis 
and experiment, the smoothing or bridging relationships described are applicable. 
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A 2:1 Ratio Inductive Voltage Divider With Less Than 0.1 
PPM Error to 1 MHz 


Cletus A. Hoer* and Walter L. Smith* 


Institute for Basic Standards, National Bureau of Standards, Boulder, Colo. 80302 
(October 25, 1966) 


\ simple ratio transformer having a 2:1 ratio of input voltage to output voltage can be made with 
a ratio error less than 0.1 ppm from 1 kHz to above 1 MHz. Applications and sources of error are 
discussed. Experimental results leading to an optimum transformer design are given. A bridge to 
measure the ratio error to 0.025 ppm from 1 kHz to 1 MHz is also described. 


Key Words: Attenuator. bridge, high frequency, inductive voltage divider, ratio transformer. 


Introduction Figure 2b shows n basic 2:1 ratio dividers connected 
to give a voltage ratio 


Techniques for winding inductive voltage dividers 
with adjustable ratios have been improved during the 
last few years to where dividers having errors less than 
0.1 ppm at frequencies around | kHz are not too dif- 
ficult to make. Unfortunately, as the frequency in- 





creases, the ratio error increases rapidly so that this 
accuracy cannot be maintained above about 10 kHz. 
At frequencies above 10 kHz the most accurate ratio 
device is an inductive voltage divider which has a 
fixed voltage ratio of 2:1: that is, one dividing an input 
voltage into two equal parts such as is shown in figure 1. 
Such dividers have recently been made with errors 
less than 0.1 ppm from 1 kHz to above 1 MHz. This 
paper summarizes the research that was done to de- 
velop this accurate high frequency divider. 


<= 
2 | 
r0209009-4 - 











FIGURE |. A 2:1 ratio inductive voltage divider. 


2. Applications 





A few of the many possible applications for this : ae 
type of divider are shown in figure 2. Figure 2a is a .> x . wee a IF re SWITCH ON (UP) 
bridge for comparing two impedances, z; and 2. is) 2 O IF i~- SWITCH OFF 


*Radio Standards Engineering Division, NBS Laboratories, Boulder, Colo. 80302. FIGURE 2. Some possible applications for 2:1 ratio dividers. 
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FIGURE 3. Typical error curve for one decade of a low-frequency 


decade voltage divider. 


Each section of this divider accurately divides its 
input voltage in half to give an attenuation of 6 dB.! 
An attenuator made with n sections will have approxi- 
mately 6ndB. 

Figure 2c shows an adjustable ratio “binary” induc- 
tive voltage divider. The principle of operation is 
the same as for low frequency decade voltage dividers 
except that each section divides its input voltage into 
only two equal parts instead of 10 equal parts. T. L. 
Zapi et al., have shown that the error curve for any 
one isolated decade of a low frequency decade voltage 
divider has an “S” shape such as shown in figure 3 |1].? 
As the frequency increases, the magnitude of the error 
increases but the curve still tends to go through zero 
at a ratio of 0.5. It is this error curve that suggested 
that 2:1 ratio dividers would be accurate at high fre- 
quencies as well as low frequencies. If n of these 2: ] 
ratio dividers are connected as shown in figure 2c the 
voltage ratio is 

n 


_< _ | | if the ith switch is up 
eS ap 


Vin 2 lo if the ith switch is down. 


The input of each section is connected either across 
the top half or across the bottom half of the preceed- 
ing section. To obtain the same resolution with a 
binary divider as with a decade divider requires that 
the binary divider have approximately 3 times as many 
sections as the decade divider has. 

Figure 2d shows three basic 2:1 ratio dividers con- 
nected in with more 2:1 ratio dividers 
connected so as to assure that the same voltage is 
applied to each of the three in series. Any number 
could be connected in series, however the combina- 
tion shown gives the very useful ratios 


series two 


a ae 
i? 1, 2, and 5. 


This type of divider could be used in a bridge such as 
shown in figure 2a to compare standards of impedance 
which are often built to have nominal values in a 1-2-5 
sequence. 


More accurately 20 log 2 


6.020600 dB 
? Figures in brackets indicate the literature references at the end of this paper. 


3. The 2:1 Ratio Divider 


The 2: 1 ratio divider is a center-tapped transformer 
that divides a voltage into two equal parts. It is con- 
venient to define the complex ratio error, p= p,+Jpi, 
as the error in the ratio of the output voltage to the 
input voltage. The output voltages e; and e2 shown in 
figure 1 then. become 

¢ 


=F (tp), 


ei 
y= — (5 — * 
> 9 ] p) 


If |p| <1, the ratio e;/e2 is approximately 
(1) 


To get an expression for the ratio error, consider the 
divider as two impedances in series which differ by a 
small amount AZ as shown in figure 4. If no current 
is drawn from the center tap, the voltage ratio e;/e2 is 


e, Zt+AZ AZ 
— = = ] + — (2) 
e» Z Z 


From (1) and (2) we get 


AZ _AR+jwAL , 
Pe 2-—SCRR +joL (9) 


where R+jwL is the open circuit input impedance, 
Zi, of the transformer. The real and imaginary com- 
ponents of p obtained from (3) are 


RAR +o*LAL 
- R2 ss w?*L? 
— wRAL = @LAR : 
oi R? T w*L? 





If Q?>1 as is often the case, where Q 


=whl/R, 
and (5) reduce to 


AL AR 
i 7, wLQ 

eee AR | AL 
” LO 


aL 








FIGURE 4. Simplified circuit of a 2:1 ratio divider. 





These equations show the importance of winding the 
transformer so that it has a large inductance. The 
value of AL/L is made small by winding the trans- 
former on a high permeability core to increase L, and 
by twisting the wires to minimize AL. A twisted cable 
is made by folding one wire in half and twisting the 
pair fairly tightly with a hand drill as shown in figure 5a. 
This cable is wound on the magnetic core and con- 
nected as shown in figure 5b. The wires are con- 
nected so that the current path from the upper tap 
to the center tap is the same as the current path from 
the center tap to the bottom tap. Thus, both windings 
link very nearly the same magnetic flux and have a 
small difference in inductance AL. Twisting the 
wires helps assure that both windings link the same 
flux. 

Figure 6 shows the variation of p with frequency 
for a transformer having 42 turns of twisted No. 26 
magnet wire wound on a ferrite torodial core. Note 
that p, is essentially constant with frequency up to 
2 MHz which is the upper frequency limit of our 
bridge. To see how high in frequency p, will remain 
constant, we need to look at the equivalent circuit of 
the divider. 

Figure 7 shows one of two exact equivalent circuits 
for an autotransformer [2] which is what the 2:1 ratio 
divider is. The Y,. is the open-circuit admittance 
measured across the input terminals and Z,,. is the 
short-circuit “leakage” impedance measured across 
either output with the input terminals short-circuited. 
It can be seen from figure 7 that the voltage ratio 
V,/V2 at the output terminals will be equal to e;/e. 
only if no current flows through Z,.. If within the 
transformer housing the admittance from A to B dif- 
fers from that from B to D by an amount AY,, and if 


Z-AY;|<1, then 


} 2 eC» 


If AY; is mainly a difference in capacitance, AC;, 


CLAMP TO KEEP LEADS 
FROM TWISTING> 


—“ ame. 
/a3tt $5 
A, \ 








/ 


B, (a) 
A, A, 


B, 
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FIGURE 6. Error curve for one of the first experimental dividers 
showing how p, and p, are relatively constant with frequency. 


V cee 
7 =1+2p+ 207L,.AC,—j2R,.wAC, (9) 


where L,. is the leakage inductance defined by Z,, 
=R,.+joL,.. A typical value of L,., such as for the 
transformer whose ratio error is shown in figure 6, is 
L..=0.1 wH. By using trimmer capacitors from A 
to B and B to D, AC; can be made less than 0.01 pF 


giving, in ppm, 


2w?L AC, < 0.1 f2 


MHz 


and 
2R,.wAC; < 0.01 Sunz 


where fz is the frequency in MHz. Since this load- 
ing error is proportional to the frequency squared, it 
will be the predominant error term at frequencies 
much above 1 MHz. The transformer whose char- 
acteristic is shown in figure 6 has trimmer capacitors 
across the output leads to adjust AC; to a minimum. 
Without adjusting AC, to a minimum the p, curve 
increases to a value of about 5 ppm at 2 MHz. 








FIGURE 5. The wires are twisted, wound on the core and connected FIGURE 7. An “exact” equivalent circuit for an autotransformer 


as shown here to make a 2:1 ratio divider. 


with lead admittances Y and Y + AY added. 
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4. A Bridge to Measure p 


A simplified diagram of the bridge used to measure 
p is shown in figure 8. The two arms of the trans- 
former are balanced against two capacitors C’, each 
having a nominal value of 1000 pF. The bridge is 
balanced with the capacitance C and the conductance 
G. The transformer is then removed from the bridge, 
rotated so that D is at the top and A is at the bottom, 
and then replaced on the bridge. The bridge is again 
balanced as before, changing C an amount AC and G 
an amount AG. Solving the bridge equations gives 


AC 


Pr=40 (10) 


— AG 


pi= lo’ (11) 


where it has been assumed that |p|<1, C” <C’, and 
G"? < wC". 

A more detailed diagram of the bridge is shown in 
figure 9a. The capacitors C and C’ are three-terminal 
capacitors having admittances to ground shown as 
dotted lines. If these admittances are unequal they 
will load the transformer causing currents to flow 
through Z,. and, thereby, change the apparent voltage 








FIGURE 8. Simplified circuit of a bridge to measure the ratio error 


of a 2:1 ratio divider. 








Dj 











FIGURE 9. A more detailed drawing of the circuit shown in figure 8. 


ratio of the transformer. To measure the true voltage 
ratio e;/e. of the transformer, a second detector, Ds, 
is placed between the center tap of the transformer 
and ground in order to balance the admittances to 
ground with Y. With both detectors nulled for both 
positions of the transformer, the bridge equations 
yield (10) and (11) as before. 

To measure the small changes in conductance, 
AG, the T-network shown in figure 9b was used for 
G. Changing C, an amount AC, gives a change in 
conductance 


AG = w#C,RAC, (12) 
and a change in capacitance which is negligible com- 
pared to AC. Putting (12) into (11) gives 


RC, 


Pi=— Fe wAC,. 


(13) 


The resolutior for p, and p; is 0.025 ppm from 1 kHz 
to 1 MHz and 0.25 ppm to about 3 MHz. A photo- 
graph of the bridge is shown in figure 10.° Bridges 
of similar design are described in references [3] 
and [4]. 


5. Optimum Transformer Design 


Using the bridge shown in figure 10, a number of 
experiments were made to determine the optimum coil 
geometry, core material, type of wire, wire sizes, and 
amount of twist of the wire. 


5.1. Coil Geometry 


Theoretically the best coil geometry would be that 
which gives a maximum inductance for a minimum 


FIGURE 10. Photograph of bridge components, showing a divider, 
with its shield removed, connected at the top of the bridge. 
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length of cable. Large inductance is needed to de- 
crease p,, and minimum cable length is needed to 
make L,. small. The leakage inductance per unit 
length of cable is given approximately by 

a 


id d 
== roses aka 1 
] 0.2 7d tan P MH/m 


(14) 


where 


P=2t+6 
t= wire insulation thickness 


d= wire diameter 


5=skin depth (< 4 


l=cable length. 


If 6 >d/2, use d/2 in place of 6. Equation (14) is 
actually the inductance per unit length of two long 
parallel thin tapes of width d spaced a distance P 
apart. By defining P and d as stated under (14), the 
results from (14) agree with measured data to within 
10 percent which is sufficiently accurate for most 
applications. Note that L,. is independent of the core 
permeability. To decrease L,. it is necessary to de- 
crease the cable length. Changing the wire diameter 
or insulation thickness has only a secondary affect 
on L,.. Measured values of L,. and R,. from 100 kHz 
to 2 MHz are shown in figure 11 for No. 24 wire. The 
leakage inductance decreases slightly due to skin 
effect. The resistance R,. increases slightly faster 
than the square root of frequency primarily due to 
skin effect. 


5.2. Core Material 


Figure 12 shows how p, is decreased by increasing 
the permeability uw. These four transformers were 
identical except for the type of core material used. 
One long twisted cable was made and cut in half. One 
half was used to wind No. 10 and the other half to 
wind No. 11. Transformers Nos. 12 and 13 were, 
likewise, made from one cable cut in half. The results 
show that p, is inversely proportional to uw. The cores 
were all very low-loss ferrite toroids with loss tangents 
of about 0.002. Cores which are very lossy do not 
give the same results. The Q for transformers with 
lossy cores is so low that the second term on the right 
of (6) may become larger than the first. The second 
term is a function of the core loss as well as the 
permeability. The core material, therefore, should 
have a high pw but not be too lossy. The core material 
can be fairly lossy if AR is small enough to make the 
second term on the right of (6) negligible. 


5.3. Type of Wire 


A number of experiments with different types of 
wire led to the conclusion that wire which has a very 
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FIGURE 11. Variation of the leakage inductance and resistance 


with frequency. 
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FIGURE 12. Ratio error versus core permeability for four trans- 
formers made identical except for core permeability. 


uniform resistivity and diameter gives the lowest AR 
values. Figure 13 shows the errors for four trans- 
formers which are identical except for the type of 
wire used. These transformers have air cores to 
eliminate the effect of magnetic core properties. 
Another reason for using air core transformers is that 
the errors are larger and easier to measure. The 
shape of the single layer solenoid type winding is 
shown in figure 14. The different types of wire used 
were formvar coated magnet wire (C), TEFLON® 
coated magnet wire (Al & A3), and copper thermo- 
couple wire (T). Note that the thermocouple wire p, 
curve is constant with frequency and quite small for an 
aircore. This means that the second term in (6) is neg- 
ligible even though the Q is fairly low. Using thermo- 
couple wire we were consistently able to wind trans- 
formers with low errors which are constant with 
frequency. 
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The thermocouple wire transformers are true “in- 
ductive” voltage dividers. The other transformers 
are more nearly “resistive” dividers than “inductive” 
dividers. The inductance of each of these four 
transformers is about 15 wH. Since p,=10 ppm for 
the thermocouple wire transformer, (6) gives AL 
=150 pH. Typical values of AL are 100 to 500 pH. 
They are constant with frequency and independent 
of wu. Expressing L as L=pN?La, (6) becomes 


__AL 
“si uN? 4a 


where N=number of turns on the core and L, is the 
geometric inductance of the core. Typical values of 
AR are 50 to 500 uO for thermocouple wire and 10 
times that for the other wire we tried. 
are all for No. 24 wire. 

It should be noted that any wire, including thermo- 
couple wire, generally gives larger errors if the wire is 
stretched. A number of transformers were made with 
wire that was stretched approximately 1 percent to 
remove slight kinks in the wire. The errors were 
significantly larger than for transformers made from 
the same spool of wire with the slight kinks left in the 
wire. It is best to purchase wire wound on large 
spools so that it is relatively free of kinks when re- 
moved from the spool. 


These figures 


5.4. Wire Size 


The experiment described above was repeated 
with No. 20 wire and also with No. 30 wire. The larger 
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FIGURE 13. Comparing results from four different types of No. 24 
wire, where Al and A3 are TEFLON ® coated magnet wire with 
single (0.0006 in) and triple (0.0015 in) insulation thickness: 
C is formvar coated magnet wire, and T is enamel coated copper 
thermocouple wire. 


wire transformers gave slightly lower p; errors than 
the corresponding No. 24 wire transformers gave. 
There was no significant decrease in the p, errors. 
The results are shown in figure 15. Transformers 
made with No. 30 wire had considerably higher errors 
than those made with No. 24 or No. 20 wire. 


5.5. Cable Twist 


To determine how much the wires should be twisted, 
a number of transformers were made identical except 
for the amount the wires were twisted. Figure l6a 
shows p, and p; plotted against the number of turns 
per unit length of the cable. The length, /; of one turn 
is defined in figure 16b. These are all air core trans- 
formers wound with No. 24 thermocouple wire. The 
dots are single layer solenoid windings whereas the 
X’s are multilayer coils wound on nylon bobbins. The 
inductance, L, of each transformer is approximately 
15 wH. As can be seen from figure 16, any twist from 
four to ten turns/inch (160 to 400 turns/meter) gives 
about the same results. Below four turns/inch the 
error increases. At about 12 turns/inch (500 turns/ 
meter) the wire breaks. In terms of the twist length, 
ly, and wire diameter, d, satisfactory results were 
obtained in the range 


0.08 < e < 0.2. 
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FiGURE 14. An air core, solenoid type winding. 
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FIGURE 15. Comparing results from four different types of No. 20 
wire where Al, A3, C, and T are as defined in figure 13. 


This equation should also be a reasonable guide for 
twisting wires other than size No. 24. 


5.6. Lead Connections 


One more important consideration in winding the 
transformers is how the leads are connected. The 
first transformers that we wound were connected as 
shown in figure 17a. The ratio errors were quite 
unstable and very sensitive to slight movement of the 
leads. To find a better way of connecting the leads, a 
number of air core transformers were made identical 
except for the lead connections. Figure 17 shows the 
three different ways the leads 
these three ways, the one shown in figure 17¢ was most 
stable and gave the smallest errors. 
in stability is probably due to the smaller mutual 
inductance and capacitance between leads. The 
lower error is probably a result of the ease with which 
the corresponding leads can be cut to the same length. 
It is important that in any arrangement of the leads 
that the leads on one side of the center tap be the same 
length as the corresponding leads on the other side. 
A difference in lead length of 0.001 in will give a AL 
of about 40 pH for No. 24 wire. 


The improvement 


5.7. Cup-Core Design 


One of the easiest ways of constructing a transformer 
with magnetic cores is to wind the cable on a bobbin 
and then clamp magnetic cup cores to the bobbin such 
as shown in figure 18. Figures 19 and 20 show the 
errors for three cup core transformers constructed like 
the one shown in figure 18. These transformers were 
made with all the optimum parameters that were found. 
The cables were made with No. 24 thermocouple wire 
twisted to six turns/inch (240 turns/meter). The cores 


were connected. Of 
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FIGURE 16. Results of air core solenoid type (dots) and bobbin type 
(X’s) transformers wound with cables twisted a different number of 
turns per inch. 
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FIGURE 17. Three different arrangements of the transformer leads. 


are high permeability, fairly low-loss ferrite cores. 
The increase in p,; and p; at frequencies above 100 kHz 
is due to the capacitance between leads not being 
equal. If these capacitances were balanced as dis- 
cussed in part 3, the errors would be less than 1 ppm 
to above 1 MHz. 

With care, the error can be made less than 0.1 ppm. 
A fourth cup core transformer was adjusted to have a 
py error less than 0.1 ppm at 100 kHz by changing the 
length of one lead slightly. Then p, was adjusted to 
less than 0.1 ppm at 1 MHz by adjusting small trimmer 
capacitors placed between the leads. As figure 21 
shows, the resulting p, error. is less than 0.1 ppm from 


30 kHz to above 1 MHz. 
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FIGURE 18. A cup-core transformer having very low ratio errors like 
that shown in figures 19 and 20. 


The cover of the shield is not shown. 





= | Hett 
3B9 -5] 145 
- 6/290 
-7|270 














0 
FREQUENCY , MHz 


FIGURE 19. Results from three cup-core transformers wound with 
Vo. 24 copper thermocouple wire twisted to six turns/inch. 


5.8. Tape-Core Transformers 


Very broad band 2:1 ratio dividers can be made with 


high permeability magnetic tape cores. Figure 22 
shows the ratio error of a divider wound on a supermal- 
loy tape core having a tape thickness of 1 mil. The 
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error for the transformers described in 
figure 19. 


FIGURE 20. Phase angle 
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FIGURE 21 
adjusted to have a minimum p, by adjusting lead length and by 
balancing the capacitance between leads. 


Ratio error of a cup core transformer which has been 
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FIGURE 22. Ratio error of a transformer wound on a supermalloy 
1 mil tape core having dimensions OD=1.150, ID=0.650 and 


H = 0.375 in. 
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core has 64 turns of No. 24 copper thermocouple wire 
twisted to about 8 turns/inch (320 turns/meter). The 
capacitance between leads has been adjusted to be 
equal. The p, error is less than 5 parts in 108 from 


1 kHz to 1 MHz. 


6. Conclusion 


The 2:1 ratio transformer is a simple, yet accurate, 
divider which is easy to construct. Using the optimum 
parameters discussed in part 5, it is easy to keep the 
ratio error less than 1 ppm from 1 kHz to 1 MHz. Air 
core transformers having no magnetic cores can be 
made with errors less than 10 ppm over a reduced 
frequency range. With extra care, p, can be made 
less than 0.1 ppm to above 1 MHz. The ratio errors 
are quite insensitive to temperature, time, shock, and 
most any change in the environment since both wind- 
ings of the transformer are influenced equally by any 
change in the environment. 


Five cup-core transformers such as described in 
part 5.7 have been combined as shown in figure 2b 
to make an accurate 30 dB attenuator at 1 MHz. 


Details of this type of attenuator together with a com- 
plete error analysis will be given in a future publication. 
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A new dual-load flow coaxial calorimeter power meter has been constructed at the National Bureau 


of Standards, Boulder Laboratories. 
of the calorimeter extends up to 4 GHz and beyond. 
an error limit of 0.38 percent. 


Designed for use as a reference standard, the frequency range 
g ) 


The power range extends from 2 W to 100 W with 


Design details, error analysis, and results of intercomparison with other standards are given. 


Key Words: Coaxial, flow calorimeter, radio frequency power. 


1. Introduction 


The increasing complexity and higher performance 
characteristics of radio frequency and other electronic 
equipment in recent years has resulted in the need 
for more accurate measurement of the rf quantities. 
In rf power measurements, for example, 1 percent 
uncertainty for measurements in industrial standards 
laboratories is often required. Formerly, uncertainties 
of 5 percent were tolerable. Because uncertainties 
are accumulated in a chain of calibrations, the uncer- 
tainties in reference standards maintained by the 
National Bureau of Standards must be less than | 
percent. 

The dual-load flow calorimeter described here is 
essentially a refinement of earlier calorimeters of this 
type [1]! and was developed to meet the need for 
greater accuracy, and extend the range of NBS refer- 
ence power standards up to 100 W. This develop- 
ment made possible the intercomparison between this 
standard and the NBS dry-load calorimeter [2] (50 mW 
to 5 W), thereby increasing the confidence in each. 


2. Theory of Operation 


The calorimetric principle has been considered the 
most accurate method for the measurement of rf power. 
This principle is based upon the first Law of Thermo- 
dynamics, or the conservation of energy. The meas- 
urement of electrical power using this principle 
depends upon the complete conversion of the electrical 


*High Frequency Electrical Standards Laboratory, National Bureau of Standards 
Laboratories, Boulder, Colo. 80301. 
! Figures in brackets indicate the literature references on page 116. 


energy as delivered by a generator into thermal energy 
in a resistive load. The heat generated in the load 
results in a temperature rise in the load and its sur- 
roundings. This temperature rise is a monotonic 
function of the input power level and may be detected, 
for example, with a thermopile located between the 
load and a reference body whose temperature is stable 
with time. 

Two types of calorimeters could be constructed to 
measure power at the levels of interest. They are the 
absolute flow calorimeter and the substitution flow 
calorimeter. (The use of dry-load calorimeters at 
these power levels, 5-100 W, is not practical because 
of the large physical size of loads required and the long 
measurement time constant.) An example of a “true” 
or “absolute” calorimetric system is shown in figure 1. 
In this system, power is measured in terms of mass, 
time, and temperature by the equation P=FcAT. In 
this equation F is the mass flow rate of the calorimeter 
fluid, ¢ is its specific heat, and AT is the equilibrium 
temperature rise of the fluid. A conservative measure 
of the uncertainty with which such a measurement can 
be made is the sum of the uncertainties with which 
F,c, and AT can be determined. The value of c is 
known very accurately for the common fluids. The 
value of AT, however, is difficult to determine accu- 
rately for low input power levels, and F can be meas- 
ured with an uncertainty no better than 0.5 percent 
at usual values of flow rate. In addition great care 
must be taken to prevent heat exchange with surround- 
ings, and a correction must be made for heating due to 
friction flow of the fluid. Therefore, the uncertainty 
in power measurement using an absolute flow calorim- 
eter is usually 1 percent or more. 

Many of the above difficulties can be overcome by 
use of the dual-load calorimeter employing d-c substi- 


11 





FLOW 
REGULATOR 


= 54 


OPEN OR CLOSED 
FLOW SYSTEM 


— 
rel 


TEMPERATURE 


i 


= 


ce 


MEASURING 
DEVICE 


FIGURE 1. 


Flow type * 


‘absolute”’ calorimeter. 





FLOW RATE 
MONITOR 











t 











FLOW RATE 
ADJUST VALVE 
| 
AUTOMATIC 
FLOW RATE 
CONTROLLER 











| GALVANOMETER 








CONSTANT 
VOLUME PUMP 


CONSTANT 
TEMPERATURE 


@ 


MANIFOLD 











7 


DIFFERENTIAL 
THERMOPILE #1 








DIFFERENTIAL 7 
THERMOPILE #2 








BALANCE 
VALVES 


RETURN 





OIL RESERVOIR 





FIGURE 2. Block 


meter 2-100 H 


tution. 


\ block diagram showing the system is given 
in 


figure 2. As the name implies, the dual-load 
calorimeter consists of two nearly identical loads con- 
nected in series or parallel to the fluid flow supply. 
Radio frequency power is applied to one load while 
d-c power is applied to the other. A differential ther- 
mopile, or other temperature sensing device, detects 
the temperature difference between the two streams on 
the downstream side of the loads when the system has 
reached thermal equilibrium. The d-c power level 
is adjusted to make the temperature difference zero, 
and the rf power is equated to the d-c power. Prior 
to the above measurement, it is necessary that the 
system be balanced by applying equal d-c power to 
each load at or near the rf power level to be measured. 
If the loads are in parallel to the fluid flow, then the 
flow rate through either one or the other may be ad- 
justed for a null at the differential thermopile output. 


diagram of dual-load flou 


, dc to 4 GHz. 


TRAP 











calorimeter power 


Thus, in the dual-load substitution calorimetric tech- 
nique, accurate knowledge of flow rate and tempera- 
ture is not required for accurate measurement results. 
In addition heat exchange with surroundings due to 
external sources is not a problem as long as it remains 
constant during the measurement period. The sub- 
stitution principle is based upon the assumption that 
heat generated by d-c or low frequency power absorbed 
by a load will have the same effect on the device used 
to sense the temperature rise of the load or a fluid sur- 
rounding it as heat generated by an equal amount of 
rf power. However, in general, equal quantities of 
rf and d-c power will not produce exactly the same 
calorimeter sensor response. Thus, an error, com- 
monly known as the rf-de substitution error, may exist. 
This error arises primarily because the rf current 
distribution in the load resistor is not identical to the 
d-c current distribution. This results in a difference 
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in the distribution of heat sources in the load in the 
two cases. In flow calorimeters, this error is mini- 
mized since the load resistor and mount structure are 
in intimate contact with a moving fluid, and the heat 
generated by the load tends to be transferred to the 
liquid regardless of the distribution of heat sources. 

The main problems associated with measuring power 
by the substitution calorimetric method are the 
relatively long measurement time constants and the 
dificulty in maintaining thermal environmental and 
flow stability over this time period. The dual-load 
flow calorimeter was designed with emphasis on 
reducing some of these problems. For example, 
by reducing the time constant of the system to a 
minimum, the need for extremely stable rf sources 
and temperature control devices was reduced. 


3. Description of Calorimeter and Operating 
Procedure 


3.1. Loads 


The loads (see fig. 3) were designed for rapid transfer 
of heat and low VSWR (< 1.03 referenced to 50 2 at 
frequencies up to 4 GHz). The load resistor is a thin 
metal film deposited by vacuum evaporation onto a 
truncated conical substrate. The load is mounted 
inside a cylindrical outer conductor. This type of 
design was proposed by D. Woods [3]. The sheet 
resistivity, p, of the film is uniform over the area of 
the cone and the resistance is given by 


R= 


27 sin 0 


In b/a (1) 


where @ is the cone semiangle, and 6/a is the ratio of 
outer to inner diameter of the coaxial line. The tem- 
perature coefficient of the resistive film is less than 
10 ppm/°C thus insuring a nearly constant value of 
R at all power levels. The characteristic impedance 
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of the section of coaxial line containing the conical 
resistor is 


Zy=2™ In b/a (2) 
2a 


where Z»= Vu/e is the wave impedance of the me- 
dium. By making p/sin 6 equal to Z», Zo=R at any 
point along the length of the load. This condition 
tends to assure the same current distribution along 
the resistor for both rf and dc. Measurements using 
a time domain reflectometer indicated that small dis- 
continuities of the order of 0.1 © did exist along the 
length of the load. Thus, it appeared probable that a 
small substitution error could exist. The upper limit 
of the error was evaluated and is discussed in the 
following section. 

The resistors are capable of absorbing up to 100 W 
of power with no significant change in their impedance 
characteristics. The physical dimensions of the loads 
were made small to reduce the measurement time. 
In order to absorb 100 W without damage to the film, 
the oil had to be circulated rapidly around and through 
the resistor body. At low levels of input power 
(2-5 W) the minimum oil flow rate was such that the 
temperature rise of the oil was approximately 5 deg. 

5 em3 ; 
15 cm /min- At 
watt 


higher levels of input power, the flow rate per watt 


This flow rate per watt was about 


7.5 cm?/min 
was reduced to about ——————— so that the tempera- 


watt 
ture rise was 10 deg centigrade. 


3.2. Temperature Detection System 


This system consists of two differential thermopiles 
(T.P. #1 and T.P. #2), of 10 junctions each, and a 
sensitive galvanometer. A_ selector switch allows 
either thermopile to be connected separately or both 
in series, to the galvanometer. The junctions of the 
thermopiles are located in the path of the moving oil 
as shown in figure 2. The oil is mixed thoroughly in 
a chamber downstream from the loads and then 
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passes around and through a plastic block in which 
the junctions of T.P. #1 are located. The close 
proximity of the thermojunctions of T.P. #1 to the 
load resistors results in rapid temperature detection 
and reduces the possibility of heat leakage before 
detection. Thermopile #2 was installed as a cross 
check against thermopile #1. Its junctions are 
mounted symmetrically inside plastic holders located 
farther downstream from T.P. #1. Any differences 
in the temperature of the oil flowing in channel A as 
compared to channel B will result in an emf generated 
by both T.P. #1 and T.P. #2. Either one or the sum 
of these emf’s is detected by the galvanometer and 
within the limits of +0.02 percent, it is possible to 
detect any balance in the calorimeter system. 


3.3. Flow System and Reservoir Temperature Control 


The flow system of the dual-load flow calorimeter 
employs precision needle valves, a flow controller 
and meter, pressure regulator, constant volume pump, 
and reservoir, all incorporated into a closed circu- 
lation system (see fig. 2). Oil pumped from the reser- 
voir passes through the pressure regulator which is 
set at 15 psi. It is then filtered and flows on to the 
flow rate meter and controller which is adjusted for 
the desired flow rate. Following the flow meter, the 
oil flow is divided into two channels, A and B, by means 
of the manifold. The flow path through each load is 
indicated by the arrows in figure 3. Oil leaving the 
loads passes through the mixing chamber, across ther- 
mopile #1, and across thermopile #2, on through the 
balance valves, to the return trap, and thence to the 
reservoir. The flow division between the two channels 
is regulated by means of the balance valves Ap—Ac 
and By-Bc. These valves were placed downstream 
from the loads because this arrangement causes a 
back-pressure which resulted in better stability and 
control and insured that the load bodies were com- 
pletely filled with oil at all times. Each valve has an 
adjustment range of 350:1 and by properly setting the 
ratio of flow between the fine and course valves, very 
small adjustments can be made in the flow division 
between channels. 

The temperature of the oil in the reservoir is con- 
trolled at 28+0.02 °C by a conventional automatic 
control circuit-cooling coil combination. This pro- 
vides a nearly infinite heat sink for the system and 
helps to reduce the measurement. time constant. 


4. Estimation of Uncertainties 


Uncertainties in measurement of power with the 
dual-load flow calorimeter were minimized by careful 
design and precision flow control and adjustment. 
The uncertainty of measurement is the difference 


between the true value and the measured value. This 
difference is usually expressed as a percentage of the 
true value. 
explanation of each are given below. 


Sources of significant error and a brief 


4.1. Flow Division Instability and Thermal Effects 


The instability in the flow division through each 
channel and the thermal drift are both reflected as a 
null shift or unbalance in the outputs of the thermopiles 
between loads A and B. As mentioned earlier, the 
effect of heat exchange with surroundings is much re- 
duced in the dual-load configuration as compared to 
the single-load “‘absolute” calorimeter. Efforts were, 
nevertheless, made to minimize and equalize heat 
exchange between the two loads and their surround- 
ings. Unequal heat exchange is not critical because 
the effect can be cancelled by proper adjustment of 
the flow rate through the individual loads during the 
d-c balancing operation. It is required, however, 
that the heat exchange be constant during the time a 
measurement is being made. After initial null (zero 
temperature difference between the oil leaving loads 
A and B) was achieved, the emf output of the thermo- 
piles was recorded over a time period much greater 
than the measurement time constant which was 
approximately 5 min. Tests at several power levels 
indicated that for periods of up to 30 min the maxi- 
mum drift was no greater than 6 wV. Since an oil 
temperature rise of 10 deg centrigrade produces a 
net thermopile output of 4000 wV, a maximum emf 
unbalance of +6 wV, corresponds to 0.15 percent 
shift in the thermal balance between the loads. 


4.2. Detection System 


The degree to which the two channels are balanced 
is a function of the overall resolution of the system. 
The sensitivity of the complete system is limited by 
the temperature fluctuations in the oil and not by the 
detection system which consists of the galvanometer 
and differential thermopile. These fluctuations limited 
the resolution to 0.02 percent. 


4.3. RF-D-C Substitution Error 


As mentioned earlier, the calorimeter loads and 
temperature detecting system were designed to mini- 
mize the rf—d-c substitution error. The loads were 
designed to provide matched terminations and thus 
insure, as nearly as possible, identical rf and d-c 
current distributions. In addition, the oil flowing 
over the resistor surfaces tends to absorb all the heat 
generated in the load regardless of the distribution 
of heat sources in the load. A small portion of heat 
may be conducted away at the points where electrical 
connection is made between the load resistor and the 
mount. Negligible heat conduction to the mount 
occurs at the input end where the center conductor 
is immersed in the moving oil stream for a distance 
of 1¥2cm. At the grounded end of the resistor, where 
it contacts the outer conductor, some heat conduction 
is possible. In order to check for a temperature dif- 
ference in this area a 10-junction thermopile was 
cemented to the outer conductor of one of the loads 
and referenced to an ice bath. “Equal” levels of d-c 
or rf power were then applied alternately to the load. 
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FIGURE 4. DLFC input line efficiency versus frequency. 


At a frequency of 4 GHz, a difference in temperature 
of 0.035 °C was observed while at 2 GHz no meas- 
urable difference was discernible. Calculations were 
made (appendix) which show that at 4.0 GHz, with a 
temperature difference of 0.035 °C, an error less than 
0.1 percent exists. 

On the basis of the above tests the rf—d-c substi- 
tution uncertainty was estimated to be 0.1 percent at 
frequencies of 1 GHz and above and 0.05 percent at 
frequencies below 1 GHz. The resolution of the meas- 
urement, due to noise and ambient temperature varia- 
tions, was approximately 0.05 percent. 


4.4. Uncertainty in Input Lines and Mismatch Losses 


In the load bodies as described in the preceding 
section, the only portion of input transmission line 
not immersed in the oil stream is a length of approxi- 
mately 3.7 em. The portion immersed in oil can be 
considered as a part of the load since the /?R loss will 
be absorbed by the oil. The loss in the remaining 
3.7 cm of input line was calculated using the following 
equation: [4] 


ar= 8.6861 Jaav 7 ere. g V we tan 5| dB (3) 


2 
1 /7e 
where w=5 Vo and w= 27f, 


€= dielectric constant, 
= permeability, 


o =conductivity of metal used for inner and outer 
conductors, 
a= inner conductor radius, 
b= outer conductor inside radius, 
f= frequency, 
!=input line length, and 
tan 6=loss tangent of dielectric. 


The results of this calculation as a function of fre- 
quency are plotted in figure 4. The uncertainty in 
the calculations of the loss factor is due to the com- 
bined uncertainty in measuring a, 6, / (length=3.7 
cm), and in the assumed values of €, w, o, and tan 6. 
The combined uncertainty of these factors could be 
as great as +35 percent at 4 GHz. This was pri- 
marily due to the uncertainty in the thickness of the 
silver plating on the inner and outer conductors of 
the line, hence an uncertainty in the value of o. For 


example, at 4 GHz the following data were obtained 
from calculations for the loss of the input line: 


Loss in input connectors...............00.008 
Loss in straight region 
Loss in taper region 

Total loss 


.0009 dB 
.0019 dB 
.0033 dB 


The loss is then 0.2+0.07 percent or an efficiency of 
99.8 +0.07 percent. Thus the power at the input con- 
nector was higher than that measured by the calo- 
rimeter by the factor 1.002 +0.0007. The rf efficiency 
factor versus frequency is shown in figure 4. At lower 
frequencies the loss in the input line decreases and 
below 10 MHz the efficiency is assumed equal to unity. 

Due to the fact that the calorimeter loads are not 
perfectly matched to Zo, part of the incident power 
will be reflected. Because a measurement of the 
incident power is usually desired, reflections cause 
the calorimeter to read low with respect to the incident 
power. The impedance of both loads was matched 
to Zo (50 ©) using a time domain reflectometer so that 
VSWR < 1.03 +0.01 at frequencies from 1.0 to 4.0 
GHz and <1.015 +0.005 at frequencies below 1.0 
GHz. A plot of VSWR versus frequency for each load 
is shown in figure 5. 

At the maximum VSWR of 1.03, 0.02 percent of 
the incident power will be reflected resulting in a 0.02 
percent error in the power measurement if the reflec- 
tion loss is neglected. For a VSWR of 1.015, the 
reflected power is only 0.01 percent of the incident 
source. Combining the maximum expected uncer- 
tainty in the correction factor for rf efficiency with the 
maximum mismatch error gives 0.09 percent above 
1 GHz and 0.04 percent below 1 GHz. 


4.5. D-C Power Measurement Error 


The d-c power substituted in the loads is calculated 
from measured values of current and voltage. The 
voltage drop across a 1 2 standard resistor in series 
with the termination is measured with a potentiometer 
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FIGURE 5. VSWR versus frequency of DLFC loads. 


to obtain the current flowing into the loads, and the 
voltage drop across the loads is measured directly 
with a digital voltmeter. The power is then found by 
the simple equation P=VI. Both V and / can be 
determined accurately to 0.01 percent giving a maxi- 
mum d-c measurement error of 0.02 percent. 

The errors discussed are believed to be the only 
significant ones in the calorimeter and measuring 
system described herein. The overall limit of error 
can be found by adding the maximum values: 


0-1 GHz 1-4 GHz 





1. Flow division stability 0.15 
2. Detection system 03 02 
3. RF-D-C substitution 05 10 
4. Line and rf reflection (mismatch 

SS Sates echo Lien Geabinsnieaiebey enn .04 09 
5. D-C measurement error .02 02 


loss 





Limit of error 


0.28% 0.38% 


5. Conclusions 


The dual-load flow calorimeter was designed and 
constructed to provide a reference standard for cw 
power measurements in the range 2 to 100 W at 
frequencies up to 4.0 GHz. A maximum uncertainty 
or error limit of 0.38 percent was achieved. This is 
a significant improvement over prior capabilities at 
these power levels. Also the frequency range of the 
NBS reference standards was extended from 1.0 
GHz up to 4.0 GHz in the power range of 5 to 100 W. 

Intercomparisons at 1, 3, and 4 GHz have been made 
between this calorimeter and the reference standard 
dry-load calorimeter which has a total estimated un- 
certainty of 0.35 percent. The intercomparisons 
included measurements at power levels of 2 W and 


4 W and the disagreement was no. greater than 0.2 
percent. 
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Equa- 
tion (3) was derived using equations in this table. 


7. Appendix 


7.1. Analysis of RF—D-C Substitution Error 


In the following analysis of the rf-d-c substitution 
error, an upper bound is arrived at which, though 
derived from approximate considerations, is felt to 
have significance. 

The total input power to the calorimeter, P;, is 
related to the power absorbed by the oil, P., and the 
power P.. absorbed elsewhere, by 

P,=P,+P-. (1A) 
Using additional subscripts, rf and d-c, to indicate 
when the input power is respectively rf and d-c (1A) 
becomes 


Prac = Poa - Prac 


Pere = Pare + Pert. (3A) 
The measuring technique used with the calorimeter 
results in 


Pasc™ Pon; (4A) 


Using (2A), (3A), and (4A), one can obtain the relation 


Prac Port | 
+—Se | ot 1]. 
Prac Bs : 


? 
Prot _ 


; (SA) 
P, dc 


Thus, if Pey=Pege, there would be no substitution 
error. However, this condition may not hold for all 
frequencies and a measure of the ratio of Pers to Pea 
is obtained from the following considerations. 

As noted in the text, there is litthe chance for heat 
losses to the external environment except possibly at 
the junction of the load resistor with the outer con- 
ductor. A thermopile connected externally between 
this area and an ice bath indicated that there was no 
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difference in the temperature for frequencies up to 
2 GHz when equal rf and d-c power were alternately 
applied to the calorimeter. However, at 4.0 GHz a 
temperature difference of 0.035 °C was measured. 
Because of the construction and nature of operation 
of the calorimeter, the principle mode of heat transfer 
to the external environment is by Conduction. Since 
the rate at which heat flows by conduction is propor- 
tional to the temperature difference between the source 
and sink, 


> 7 —— 
i 


T 
amb 


: Son gn (6A) 
Prac lj 1 amb 


where 7, and T,. are the respective temperatures at 
the junction of the load resistor and outer conductor 
when rf and d-c power are alternately applied to the 


247-602 O-67—3 


calorimeter. The ambient temperature is noted by 
From experimental data at 4 GHz, 
Pert 


———— = ] 023. 
Prac 


amb: 
(7A) 


An upper bound for the ratio of Pegc/Prac, obtained from 
using the dual-load calorimeter as an absolute flow 
calorimeter, is 


Fede — 9,03. 


(8A) 

Prac 
Substituting (7A) and (8A) into (5A), the upper limit 
for the substitution error is 0.07 percent which was 
increased to 0.1 percent because of the approximations 
used. 


(Paper 71C2-250) 
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The results from a study of the effects of aging and humidity change on certain polymer specimens 


indicate it is possible to establish stable dielectric reference specimens. 


Specimens of polyethylene, 


polystyrene, polycarbonate, poly(tetrafluoroethylene) (PTFE), fluorinated ethylene-propylene polymer 


(FEP), and poly(1,4-cyclohexylenedimethylene terephthalate) were used in the investigation. 


Slow 


changes in the dielectric properties. were observed on some specimens, over a period as long as three 


years. 


Very long timed humidity runs indicate PTFE, and FEP to only a slightly less degree, are best 
suited for dielectric reference specimen preparation. 


Key Words: Dielectric constant, dissipation factor, frequency, fluorinated ethylene-propylene, 
humidity, polycarbonate, polyethylene, poly(1,4-cyclohexylenedimethylene tereph- 
thalate), polystyrene, poly(tetrafluoroethylene), time, reference specimens. 


1. Introduction 


A long term study of certain polymer specimens 
under consideration as dielectric reference specimens 
has been conducted. Such reference specimens may 
be used to compare the measurement accuracy of 
various other laboratories, so these specimens must 
have accurately determined and stable dielectric 
properties. Specimens from several polymers were 
found to exhibit sufficient stability to permit their use. 
The materials investigated were polyethylene (PE), 
polystyrene (PS), polycarbonate (PC), poly(tetrafluoro- 
ethylene) (PTFE), fluorinated  ethylene-propylene 
polymer (FEP), and poly(1,4-cyclohexylenedimethylene 
terephthalate) (T-16). All data were taken on flat 
disk specimens, and the conditions under which the 
measurements were made were carefully controlled. 
The humidity of the surrounding atmosphere as well 
as the specimen temperature was controlled because 


the properties of most specimens were functions of 


humidity. 


2. Measuring Equipment 


A low-voltage conjugate Schering bridge similar to 
one described in a paper by Scott and Harris. [1]! was 
used to measure dielectric constant and dissipation 
factor from 10? Hz to 10° Hz. The cell used with the 
Schering bridge was of the guard-ring micrometer- 
electrode type [2, 3, 4] and was mounted in a chamber 


*Present address: 4220 Franklin Street, Kensington, Md. 20795. 
1 Figures in brackets indicate the literature references at the end of this paper. 


whose temperature was controlled by pumping liquid 
from a temperature bath through coils around the 
chamber, and whose humidity was maintained by 
exposing inside the box either P.O; for < 1.5 percent 
relative humidity (rh) or a saturated Mg(NOs)2 solution 
for about 52 percent relative humidity. The humidity 
in the chamber was measured with an electric 
hygrometer. 

- In the frequency range from 1 MHz to 10 MHz a 
Q-meter was used for dielectric constant and dissipa- 
tion factor measurements. The cell used was a two- 
terminal micrometer-electrode cell mounted in a box 
whose temperature and humidity were controlled in a 
manner similar to that used for the three-terminal cell. 
A two-terminal Schering bridge was used to coordinate 
the Q-meter measurements with the more accurate 
three-terminal data. 


3. Specimens 


3.1. Contact Electrodes 


Many of the measurements were made without con- 
tact electrodes using the air-gap technique [5] with 
8.5 cm diam disks from 2 mm to 5 mm thick. Disks 
with 3.8 cm diam were prepared with contact elec- 
trodes so that specimens might be measured at higher 
frequency with the Q-meter and used to calibrate 
two-terminal equipment. Heavy gold _ electrodes, 
about 1500 A thick, were evaporated on the faces of 
these similar specimens to form the contact electrodes. 
Care was always taken to see that the resistance across 
the electrode surface was never more than a few tenths - 
of an ohm. 
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3.2. Surfacing Specimens 


At the beginning of this study a modified Schiefer 
abrasion machine was used to surface the specimens 
[6]. Even though this method produced specimens 
with thicknesses uniform to within 6 w, we later used a 
technique which produced more uniform thicknesses. 
The second method made use of a vacuum chuck and a 
fly cutter on a milling machine. The vacuum chuck 
consisted of a metal vacuum chamber with holes drilled 
in one side for holding the specimen in place and could 
be positioned on a milling machine for successive 
passes made with the cutter. The specimen was 
turned over or rotated with each pass until the faces 
became quite parallel. With careful preparation of 
the cutting tool and proper adjustment of the cutting 
rate, this technique could produce specimens with 
thicknesses uniform to within 3 yw over an 8.5 cm diam. 


3.3. Measuring Thicknesses 


The specimen thicknesses were measured with an 
indicating micrometer having a reference point set by 
gage blocks. Measurements were made at points 
around concentric circles as determined by a pattern. 
These numbers were averaged to determine the meas- 
ured thickness. “omparisons indicate agreement to 
0.05 percent between the micrometer thickness and 
that determined with the two-fluid electrical technique 
[7]. 

3.4. Storage 


Between measurements and for initial conditioning, 
specimens were stored in glass vessels with humidity 
controlled in the same manner as in the measuring 
chamber. 


3.5. General Properties 


The suppliers indicated the following general proper- 


ties for the various materials. The polyethylene was 
the high molecular weight, low pressure type. It had 
a linear structure with no detectable branching and a 
crystallinity of from 90 to 95 percent. The T-16 poly- 
ester was a linear, thermoplastic, crystallizable poly- 
mer about 36 repeating units long. It had a crystal 
melting point of around 265 °C and a glass transition 
at about 87 °C. The polycarbonate material con- 
tained no plasticizers and had a molecular weight 
ranging around 30,000 to 35,000. The polystyrene 
material was reported to be of high purity. PTFE and 
FEP specimens were from laboratory grade stock. 

The polyethylene specimens and certain of the poly- 
styrene specimens were compression molded in our 
laboratory. The polycarbonate specimens were im- 
pact molded, the T-16 specimens were injection 
molded, and the FEP and PTFE specimens were com- 
pression molded, all by the supplier. 


3.6. Density 


The density of most specimens was measured in 
order to characterize the materials better and to pro- 
vide a rough measure of the uniformity of specimens 
from the same material. Two methods were used for 


determining density. In the first method the volume 
of the specimen was determined by liquid displace- 
ment. The second method was based on dimensional 
volume measurements. Since our specimens were 
very carefully machined disks, they were, in fact, good 
approximations of right circular cylinders. By meas- 
uring the diameter with a traveling microscope and 
using the thickness, we could calculate a dimensionally 
determined volume. The same mass data was used 
with each volume to calculate two different densities. 
Based on repeated measurements and calibration of 
the traveling microscope, we estimate the dimensional 
volume measurement to be accurate within about 0.1 
percent; however, the liquid displacement densities 
were consistently larger than the dimensional densities 
by 0.1 percent to 0.5 percent, even though care was 
taken to reduce errors caused by insufficient surface 
wetting of the specimen and by surface tension acting 
on the suspension wire. Since the dimensional meas- 
urements also produced much more reasonable rela- 
tionships between the dielectric constant and density, 
it was concluded that our liquid displacement measure- 
ments were less accurate, so all the density data pre- 
sented were determined from mass and dimensional 
measurements. 


4. Measurement Procedure 


The following procedure was used to obtain the 
dielectric properties of individual specimens as func- 
tions of time after humidity change. Each specimen 
was stored in a chamber with the desired relative 
humidity until it reached equilibrium, i.e., until its 
dielectric properties no longer changed with time. A 
timed humidity run was then begun by transferring 
the specimen to another chamber at the new humidity. 
Periodic measurements were made with careful note 
being taken of the elapsed time since the humidity 
change. All timed humidity run data as well as fre- 
quency sweep data were taken at 23 °C. 

It was observed for some of the electrodeless speci- 
mens that slightly different values of dielectric constant 
and dissipation factor were obtained when the speci- 
men was turned over. The amount of this effect 
varied from specimen to specimen, and for some disks 
only the dielectric constant differed significantly. 
Differences as large as 0.2 percent in the dielectric 
constant were observed. This behavior was caused 
by the specimens being slightly curved instead of flat 
although they were of uniform thickness. In some 
cases it was possible to partially flatten out the warped 
specimens, and this reduced the difference. When- 
ever possible, the same face was kept turned up 
throughout the timed humidity runs; however, for some 
of the older data, taken before the effect was dis- 
covered, this was not done. 


5. Results 
5.1. Timed Humidity Runs 


Figures 1 through 7 show timed humidity runs meas- 
ured at 1000 Hz for certain specimens of the various 
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materials. For all of these figures, time from humidity 
change is plotted on the horizontal axis using a logarith- 
mic scale, and those points marked with arrow heads 
on the left-hand margin of some of the plots indicate 
zero time starting values, which are measurements 
made at the storage humidity shortly before the begin- 
ning of the timed humidity run. Even though there 
are usually curves from only one or two specimens for 
a given material, the trends are generally representa- 
tive for several of our specimens. All the dielectric 
properties given here were measured using the air- 
gap technique except when contact electrodes are 
indicated. 

Figure 1 shows a plot for a T—-16 specimen going 
from <1.5 percent rh to 52 percent rh after being 
exposed to the low humidity for 455 days. This 
material does not show long term stability as indicated 
by the negative slope of the curves after about 100 days. 
We suggest that this aging process may be due to very 
slow crystallization of this material. We offer this as 
a possible explanation, even though the specimens 
were kept well below the glass transition temperature, 
because the material was not completely amorphous 
and the trend is indeed in the direction of higher 
crystallinity. Other specimens of this material in 
both humidities also showed this downward trend 
with time. 

Figure 2 is a plot for two contact electrode speci- 
mens of polystyrene going from 52 percent rh to 
<1.5 percent rh. The specimen for curve A had been 
in 52 percent rh for 1600 days, and the specimen for 
curve B had been in 52 percent rh for 74 days. 

In figure 3 are curves for two specimens of poly- 
carbonate. The specimen producing curve A was in 
52 percent rh for 992 days before the humidity change 
and likewise the specimen for curve B was in < 1.5 per- 
cent rh for 1096 days before the change. These curves 
show the reversibility of the process, but the long time 
data seem to indicate an aging process again. Figure 
4 shows curves for a polycarbonate specimen differ- 
ently prepared by the manufacturer; this specimen 
had been in 52 percent rh for 174 days. Aside from 
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FIGURE 2. Timed humidity runs for two contact electrode specimens 


of PS going from 52 percent rh to < 1.5 percent rh. 
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FIGURE 3. Timed humidity runs for two specimens of PC. 


Specimen A going from 52 percent rh to 1.5 percent rh. and specimen B going from 


1.5 percent rh to 52 percent rh. 
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FIGURE 4. Timed humidity run for a specimen of differently pre- 


pared PC going from 52 percent rh to < 1.5 percent rh. 


a lower dielectric constant and dissipation factor, the 
behavior indicated is about the same. 

The data in figure 5 are for a specimen of FEP which 
was in < 1.5 percent rh for.448 days before the change; 
these curves show a slight upward trend going from 
< 1.5 percent rh to 52 percent rh. However, within the 
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FIGURE 5. Timed humidity run for a specimen of FEP going from 


< 1.5 percent rh to 52 percent rh. 
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FIGURE 6. Timed humidity run for a specimen of PTFE going from 


< 1.5 percent rh to 52 percent rh. 
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FIGURE 7. Timed humidity runs for three specimens of PE. 


Specimen A going from < 1.5 percent rh to 52 percent rh, and specimens B and C going 
from 52 percent rh to < 1.5 percent rh. 


limit of our estimated accuracy, no change is indicated 
by the curves for the PTFE specimen in figure 6. 
This specimen was in < 1.5 percent rh for 453 days 
before the change to 52 percent rh. In fact, no change 
was observed for any specimen of PTFE. 

Figure 7 contains curves for specimens of poly- 
ethylene. Curve A shows the very slow rise in dissi- 
pation factor going from < 1.5 percent rh to 52 percent 
rh for a contact electrode specimen which had been 
in < 1.5 percent rh for 199 days. No change in dielec- 
tric constant has been observed for any polyethylene 
specimen. The specimen producing curve B was in 52 
percent rh for 700 days before the change while the 
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FIGURE 8. An expansion of timed humidity runs for three specimens 
of PC going from 52 percent rh to < 1.5 percent rh. 


contact electrode specimen associated with the curve 
C was in 52 percent rh for 123 days. Taken together 
these data would seem to indicate that, even though 
water is sorbed into polyethylene slowly, some of the 
water is released very quickly. 

In figure 2 there are peaks in the PS dissipation 
factor curves after one day. Figure 8 shows similar 
peaks for three specimens of PC on an expanded time 
scale, but for this material the peak is in the region 
of 10 to 30 days. This strange behavior was observed 
only on PS and PC specimens. Moreover, it has been 
observed for specimens measured with contact elec- 
trodes as well as with the air-gap method. No 
explanation has been found for this unusual data. 


5.2. Frequency Sweeps 


Figure 9 contains the dielectric constant and dissipa- 
tion factor plotted against frequency for a contact 
electrode specimen of PC; these data were measured 
at <1.5 percent rh equilibrium. Figure 10 likewise 
shows the dissipation factor and dielectric constant 
for contact electrode specimens of PS at both < 1.5 
percent rh and 52 percent rh equilibrium. Figure 11 
contains similar plots for two T—16 specimens. 

From the other materials —PE, PTFE, and FEP —no 
specimens showed a significant change in dielectric 
constant over our frequency range. Dissipation factor 
data for FEP specimens at both the low and high 
humidity appear in curves A and B in figure 12. The 
dissipation factors for PE specimens at < 1.5 percent 
rh measured less than the limit of our accuracy so 
were not plotted, but 52 percent rh data appear as 
curve C in figure 12. 

Results differed slightly from specimen to specimen 
for most measurements; this was especially true at 
52 percent rh. We were ultimately concerned more 
with the careful recording of the behavior of a particu- 
lar specimen or group of specimens than with the 
general material, so none of the data presented here 
should be taken to represent the whole material in 
a detailed way. 
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FIGURE 9. Dielectric constant and dissipation factor versus fre- 
quency for a contact electrode specimen of PC at < 1.5 percent rh 
equilibrium. 
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FIGURE 10. Dielectric constant and dissipation factor versus fre- 
quency for two contact electrode specimens of PS. 
) 


Specimen A is in 52 


percent rh equilibrium and specimen B is in < 
equilibrium. 


1.5 percent rh 


5.3. Density and Water Absorption 


A linear relationship between the dielectric constant 
and density is a good approximation over a restricted 
density range for polyethylene. The common form 
for this relation is the Clausius-Mossotti equation, 


nl 


E 2 


where p is the density, €’ is the dielectric constant, and 
K is a constant for the material which depends on the 
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FIGURE 11. Dielectric constant and dissipation factor versus fre- 


quency for two specimens of T—16. 


Specimen A is in 52 percent rh equilibrium and specimen B is in < 


1.5 percent rh 
equilibrium. 
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FIGURE 12. Dissipation factor versus frequency for specimens of 


FEP and PE. 


Specimen A is FEP at < 1.5 percent rh equilibrium; specimen B is FEP at 52 percent rh 


molecular weight and the total polarizability per mole- equilibrium; and specimen C is PE at 52 percent rh equilibrium. 
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FIGURE 13. Dielectric constant versus density for specimens of 


polyeth ylene. 


cule. The dielectric constant of six PE specimens is 
plotted versus density in figure 13. A straight line 
was fitted to the data by a method, useful when the 
‘errors in the variables are unknown, as given by 
Bartlett [8], and this line is also shown. These dif- 
ferences in density were probably caused by the 
slightly different thermal and mechanical treatment 
received during molding and machining. 


Some effort was made to determine the amount of 


water absorbed by certain specimens by taking the 
difference between the equilibrium weight at < 1.5 per- 
cent rh and weight at room humidity which was gen- 
erally 20 percent rh to 40 percent rh and never above 
50 percent rh. A difference in weight of 0.03 percent 
was obtained from the average change of two PS speci- 
mens. This roughly agrees with the value observed 
by von Hippel and Wesson [9] of 0.048 percent for 
60 percent rh at 25°C. Likewise, for PC we observed 
a 0.13 percent change in weight which compares with 
the 0.15 percent observed by Thompson and Goldblum 
[10] for 50 percent rh. These can not be exact com- 
parisons since the materials were, at best, only similar 
in composition, but they do show order of magnitude 
agreement. We were unable to detect a difference in 
weight for FEP, PTFE or PE even though the dissipa- 
tion factor of PE showed a quite definite change going 
to 52 percent. The amount of water absorbed by PE 
has been shown to be quite small (on the order of 0.03 
or 0.06 mg/cm*) [11], 12], and beyond the sensitivity 
of our weighing system. 


6. Estimated Error Limits 


From the use of this low-voltage conjugate Schering 
bridge over a long period of time, it is estimated that 


repeated measurements of capacitance over the fre- 
quency range 100 Hz to 100 kHz could be expected 
to agree to within +0.001 pF for measured capacitances 
above 10 pF and 0.0001 pF for 10 pF and below. 
Based on this estimate and comparisons with the two- 
fluid technique as well as repeated determinations, we 
estimate our dielectric constant measurements using 
this bridge and the air-gap technique to be accurate 
within + 0.1 percent. 

Previous papers by Scott and Harris [2, 13] discuss 
some of the important problems involved in the meas- 
urement of small dissipation factors. Significant 
errors can arise if both the high and low sides of the 
circuit have some finite common impedance to ground. 
In general, this is the case with a three-terminal cell 
or capacitor. It is possible, in principle, to completely 
analyze the circuits involved for a given cell and bridge 
to allow the most accurate dissipation factor deter- 
mination: however, that study has not yet been com- 
pleted for our equipment. Heating by the measuring 
field can also be a problem with higher loss specimens, 
since the properties may change rapidly with tempera- 
ture and frequency, so we have included large per- 
centage limits with the estimated accuracies for the 
dissipation factor. 

From the use of this Schering bridge over a long 
period of time, in the same manner as capacitance, it 
is estimated that repeated measurements of the dis- 
sipation factor of a stable specimen could be expected 
to agree within 2 X 10~-®, and based on this and com- 
parisons with carefully constructed low loss capacitors, 
we estimate the absolute and relative dissipation factor 
measurements using this bridge to be accurate within 
the greater of the two limits given in their respective 
columns in table 1. 

The dissipation factor measurements made using the 
Q-meter with an extra-high-Q coil at 1 MHz are 
estimated to be accurate within +0.0002. At 10 
MHz, this accuracy was reduced somewhat due to the 
uncertainty of corrections for lead inductance. The 
dielectric constant measurements made with this 
Q-meter are believed to be accurate within +0.2 
percent. These estimates were derived from measure- 
ments made on specimens whose properties at 1 MHz 
and 10 MHz were believed to be essentially the same 
as at lower frequencies, where they could be deter- 
mined more accurately. 


TABLE 1. Estimated error limits for absolute and relative dissipation 
factor measurements made with the low-voltage conjugate Schering 
bridge are given by greater of the two values 


Frequency Absolute error limits Relative error limits 


102 + 30 x 10-* or + 2% 
10° +20 x 10-* or 
10* +20 x 10-* or 
10° +40 x 10-* or 


+20 x 10-* or +2% 
+10 X 10-* or +1% 
+10 10-* or +1% 
+30 x 10-* or +3% 


7. Discussion 


The PC and T-16 timed humidity runs show large 


rapid changes as well as the long term drift. If these 
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materials were to be used for making dielectric refer- 
ence specimens, it would be necessary to carefully 
control the humidity, and obviously their usefulness 
would be limited at best. The PC data does give a 
clear indication of the reversibility of this process for 
the same or a similar specimen, and therefore, if a 
copy of the timed humidity run were sent along with a 
specimen, a somewhat more accurate characterization 
could be obtained under a large variety of conditions. 

For the PS specimens, the change was not nearly 
so great; however, it would still be necessary to take 
the humidity into account to obtain an accuracy better 
than about | percent for dielectric constant and about 
1 X 10-4 for dissipation factor. 

The PE data show a very rapid initial change in the 
dissipation factor when the relative humidity is 
changed from 52 percent rh to < 1.5 percent rh, but 
quite a slow drift upward on the return to 52 percent rh. 
This behavior coupled with the lack of a change in the 
dielectric constant means that it would be possible 
to predict a specimen’s properties quite well if it were 
stored in low humidity before and after a measurement 
even if the measurement were made at a_ higher 
humidity. 

PTFE, and FEP to possibly a slightly less degree, are 
materials well suited for dielectric reference specimen 
preparation. Although the FEP dielectric constant 
and dissipation factor data seem to show an upward 
trend going from <1.5 percent rh to 52 percent rh, 
lines fitted by the least squares method versus the 
logarithm of time after humidity change do not have 
slopes statistically different from zero. For PTFE, no 


change in dielectric constant or dissipation factor is 


indicated within our error limits for the 100 Hz to 
10 MHz frequency range or for a change in humidity. 
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The Sensitivity of the Dicke Radiometer* 
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The literature is reviewed concerning the sensitivity of the Dicke radiometer, excluding gain fluc- 
tuations. Discrepancies are pointed out and a new derivation of sensitivity using a Fourier transform 
method is used to resolve these discrepancies and to extend the results to radiometers with lossy 
switches. 

Experimentally it is shown that radiometers using a half-wave square-law, linear-law, intermediate- 
law, or envelope detectors all have a sensitivity equal to the theoretical full-wave square-law detector 
(within the +20 percent uncertainty of the experiment). 


Key Words: Dicke, microwave, radiometer, sensitivity, switching. 


1. Introduction 


At the National Bureau of Standards a switching radiometer is used to compare microwave 
noise sources with a standard hot load [1],! and a sum-and-difference-correlation radiometer is 
being developed to compare microwave cryogenic noise sources which uses a switching radiometer 
as a component [2}. 

Although there is extensive literature discussing the Dicke radiometer, existing publications 
are dificult to use. This is partly due to a series of errors and corrections and partly due to 
divergent assumptions, notation, and forms. The literature as a whole has unresolved differences 
between various statements of the sensitivity for the Dicke radiometer. 

The purpose of this paper is to review the literature (sec. 2), derive the sensitivity of the 
Dicke radiometer in a form which includes switch noise (sec. 3), and to compare the derived 
sensitivity with experiment (sec. 4). The sensitivity is rederived, rather than making a further 
correction, for two reasons. First, a statement of sensitivity is needed in a form general enough 
to include the prior literature. This facilitates identifying discrepancies. Second, a basically 
different derivation helps to avoid prior errors and to increase confidence where agreements are 
obtained. 


2. Review of the Literature 


To facilitate comparison between the various authors in the literature, a brief discussion of a 
model switching radiometer is presented. The literature is compared in terms of this model. 
The symbols of this paper may be converted to the notation of the other authors by means of 
appendix A. Also contained in appendix A are equivalent statements of predetector bandwidth 
B, and postdetector integration constant T. 


*Sponsored by Advanced Research Projects Agency, Project DEFENDER; ARPA Order No. 515. 
**Microwave Circuit Standards, National Bureau of Standards Laboratories, Boulder, Colo. 80302. 
! Figures in brackets indicate the literature references at the end of this paper. 
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2.1. Basic Switching Radiometer 


The basic switching radiometer is shown in figure 1. Essentially the power from a noise 
source at the effective temperature 7, is compared with a noise source at effective temperature T». 
It is usually advantageous to operate the radiometer with a known and variable JT, which is made 
equal to 7;.. The comparison of T; to T2 is effected by switching the input of a receiver from 
T, to Tz and examining the output of the receiver for a signal correlated with the switching fre- 
quency. Uncorrelated terms (for example, the internal noise of the receiver) do not contribute 
to the de output of the correlator. 

The noise voltage y;(t) from the noise source at temperature 7; is connected to the receiver 
through a variable coupler of voltage transmission coefficient a;(t). Similarly y(t) from noise 
source 7. is fed through a variable coupler with voltage transmission coefficient a2(t). The internal 
noise of the amplifier referred to the input of the amplifier has effective temperature 7,. The 
signal y,(t), the input to an ideal noiseless amplifier, is 


ys(t)=ay(t) yi (t) + ao(t) y2(t) + yn(t). (1) 


where y,(t) is the voltage needed to represent the amplifier noise. Ideally the transmission 
coefhicients a,(t) and a2(t) vary from zero to one periodically with period vo!. Typically ai(t) and 
a2(t) are 180° out of phase. The transmission coefficient a;(t) is referred to as the voltage modula- 
tion. The power transmission coefficient p(t) = a(t) is referred to as the power modulation. 
From the modulation a;(t) is generated a periodic reference signal c(t) which then goes to the multi- 
plier. The reference signal c(t) is referred to as the correlation and usually is sinusoidal or square 
wave. 
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2.2. General Comments 


Eleven papers (as noted in table 1) are here reviewed which contain statements concerning 
the sensitivity of a Dicke radiometer. There is an underlying agreement among eight of these 
and the remaining three retain many of the essential features. In order to recognize this agree- 
ment, various definitions must be related (see appendix A) and account taken of the modulation 
and correlation assumed. 

There is an important difference among the various agreeing papers. The nature of this 
difference is essential to the general understanding of the literature, so it will be discussed before 
commenting on individual papers. 

The radiometer described by Dicke [3] modulates the noise source by means of a variable 
attenuator which has nearly the same temperature as the noise source of interest. The papers 
of Goldstein |4, 5], Strom [6, 7], Strum [8], Knight [9], and Johnson [10], collectively referred to as 
GSSKJ, calculate the sensitivity of a somewhat different radiometer. The GSSKJ version of the 
Dicke radiometer modulates the noise source of interest by an ideal noiseless device, which could 
be a variable attenuator at the absolute zero of temperature. This corresponds to T, = 0 °K in 
the model introduced in section 2.1. The Dicke radiometer was originally a null device but this 
concept was extended by Goldstein to include an off-null operation. This introduces many prac- 
tical differences. For example, a radiometer operated off-null is not immune from receiver gain 
instabilities. The radiometer described by GSSKJ is a special case of the papers of Kelly, Lyons 
and Root [11], 12], Colvin [13], and Tiuri [14] collectively referred to as KLRCT. However, GSSKJ 
and KLRCT differ in detail. For example, GSSKJ’s results indicate somewhat greater fluctua- 
tions in the radiometer output due to the input noise source (7;) than the corresponding KLRCT 
results. 

The origin of the discrepancy lies in the different technique of obtaining the output of the 
square law detector relative to the input. GSSKJ use an autocorrelation method and handle the 
resulting squared terms by means of a theorem usually attributed to Rice [15]. Rice’s theorem 
for square law devices is rigorous for stationary Gaussian processes. KLRCT essentially apply 
Rice’s theorem directly, neglecting the usually small nonstationary aspect of the signal at the detec- 
tor input. The errors introduced by this simplifying assumption go to zero as the noise sources 
being compared approach the same temperature (7;=T2), and the receiver is switched between 
the two sources in such a way that the power to the receiver is constant throughout the switching 
cycle. 

The calculations in this paper are as precise as the GSSKJ calculations. Thus detailed 
agreement with KLRCT is expected only when 7, =72. 

In the last column of table 1, which appears later in this paper, prior calculations of the radi- 
ometer sensitivity are indicated by modulation and correlation type. Those calculations which 
differ with this paper are marked with an asterisk. Although Goldstein [4, 5], Strum [8], and 
Knight [9] are marked with an asterisk, they are considered to be in the group of eight having under- 
lying agreement. More will be said of this in the detailed comments which follow. 


2.3. Specific Comments 


In 1946, Dicke [3] suggested modulating a noise signal under study as a means of distinguishing 


the signal of interest from the noise originating in the amplifying system and also to reduce the 
fluctuations in the radiometer output due to receiver gain instabilities. He derived the sensitivity 
of a radiometer using a statistical approach for square-wave modulation and sinusoidal correlation. 
In the analysis Dicke assumed that the detector was square-law, then he generalized to the linear 
detector. Selove [16] corrects an error in the generalization to the linear detector and Strom 
[6, pp 37 and 68] challenges an overall factor. For a square-wave modulation and sinusoidal 
correlation, Dicke obtained a sensitivity which is a factor of 2 greater than the result to be obtained 
in this paper and thus comparable with the accepted sensitivity of an ideal (no gain fluctuations) 
total power radiometer [11], 12, 13, or 14]. Because the switching radiometer cannot have com- 
parable sensitivity, Dicke’s result must be in error. 
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In 1955, Goldstein [4] analyzed the switching radiometer using sinusoidal-voltage modulation 
(or switching), sinusoidal correlation and square-law detection making use of the autocorrelation 
function and assuming the noise signal had a Gaussian distribution. However, as noted above, 
he assumed a noiseless modulator which corresponds to assuming T; =0 °K. Goldstein’s result 
as corrected by himself [5] and Johnson [10] agrees with a special case of the analysis to be 
presented. 

Also in 1955 Bunkin and Karlovy [17] analyzed the switching radiometer using the autocorrela- 
tion method. Their analysis assumes a Gaussian noise distribution, square-wave modulation and 
sinusoidal correlation. In addition, the effects of amplifier gain fluctuations are considered. 
This paper depends to some extent on other papers in the Russian literature which are not readily 
available. It was not established for certain if the paper agrees in all details, but no obvious 
discrepancies occur. 

Strom |6, 7] shows that the sensitivity of the switching radiometer is unaltered if a biased diode 
is used as the detector instead of a square-law detector. He obtains a sensitivity that is less by a 
factor V2 than the expected value because the out-of-phase rejection of noise by the correlator 
was not taken into account. 

Strum [8] discusses the effect of fluctuations in various components of a radiometer system 
and discusses various types of linear detectors. He comments on the use of an IF amplifier in an 
analysis. Strum defines a generalized bandwidth which reduces to the usual one for rectangular 
bandpass but is not equivalent in general to that used by Dicke: Kelly, Lyons, and Root: Ward and 
Richardson; Tiuri; and in this paper. 

Kelly, Lyons, and Root [11], 12] analyzed a general comparison radiometer (7; and 72 arbitrary) 
for arbitrary modulation and correlation. They treated gain fluctuations and optimized the modu- 
lation wave form depending on the rms amplitude of the gain fluctuations. They also presented 
a proof that a square-law detector is superior to any other detector law for Gaussian noise (ignoring 
gain fluctuation of the amplifier). Except for the differences off-null already noted, there is no 
disagreement with this paper. 

Colvin [13] discusses the switching radiometer under rather general conditions as part of an 
extensive survey of various radiometers. Again, except for the off-null differences discussed, there 
is no disagreement with this paper. 

Knight [9] considers a radiometer with asymmetrical rectangular modulation and correlation. 
The calculation is of the Goldstein type. In his paper, eqs (3) and (4) lack the factor 27 on the 
right-hand side which occurred in his prior equations. Knight's corrected equations are a special 
case of the analysis to be presented. However, his extension to the balanced radiometer operation 
did not take into account the fluctuation of the output due to the noise sources. 


Ward and Richardson [18] consider a power sinusoidal modulated radiometer with IF ampli- 


fiers of different image and signal responses. The sensitivity obtained is that expected for a total 
power radiometer, and for the same reason as noted for Dicke’s paper, this must be in error. 
Johnson [10] corrected Goldstein’s results and extended the calculation to include various 
combinations of sinusoidal and square-wave modulations and correlations. However, in Johnson’s 
extension to square-wave modulation his (1/4) o?03, term should be replaced by o203.2 
Tiuri [14] states the radiometer sensitivity for several modulations and correlations. The 
results are the expected ones within the approximations already discussed. 


3. Analysis 


The sensitivity of the switching radiometer depicted in figure 1 is determined by calculating 
the ratio of the de output of the radiometer y-2 and the rms fluctuations at the output, op». The 
time dependent noise voltages will be denoted by y(t) with a subscript denoting the point in the 
circuit shown in the circles in figure 1. An IF amplifier is not used in figure | but is treated in 
appendix C. 


? This has been acknowledged by Johnson in a private communication. 





3.1. The Fourier Transform 


The Fourier transform preserves required phase information. Because the transform of a 
random function does not converge, all the noise voltages y(t) are defined to be zero except for 
—T/2<t<T/2. The truncated function of y(t) has identical properties with y(t) in the limit 
T—>«. Thus by y(t) it will be understood the truncated function of y(t). The Fourier transform 
of y(t) is denoted Y(v), again with a subscript to denote point in the circuit. 

The following form is used for the Fourier transform pair: 


Yon= [7 emmy 


v(t) =| ; e?ttY (py) dv. (3) 


The Fourier transform or inverse transform of a product (in either ¢ or v space) is the convolution 
of the transforms of the two factors (convolution theorem). With the choice in (2) and (3), the 
convolutions (denoted by ¥&) are defined by 


x 


veltyaeyste= | yzr(t—t’)y2(t')dt’, 


x 


Yow eY.00)= | Y,(v—v’')Y.(v')dv'. 


xs 


3.2. The Power Spectral Density and the Wiener-Khintchine Theorem 


The power spectral density W(v) is related to the Fourier transform of y(t) in the following 
manner [19], 


Lim Y( v)Y*(v 


W(v)=7, T (6) 


where the asterisk denotes the complex conjugate. It can be shown that W(v)dv is the power 
dissipated in a 2- resistor due to the spectral components between v and v+ dv in the voltage y(t). 
The Wiener-Khintchine Theorem is stated [20] 
W (v) =| R(t)e-227"'dt, 
so that 


Rit) =| W (v)ei27""dv, 


where R(t) is the autocorrelation function defined by 


Rit) = Lim = yit')yit’ + t)dt’. 
T J_qy 


Using (8) and (9) for t=0, 
v= |" W (v)dv, 


where the bar indicates the time average in the interval — 7/2 to + 7/2. 
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The analysis is ultimately stated in terms of measurable properties of the input noise signals 


(21, 22]: 


] 
W i(v) oe kT,, 


W(v) =; kT», 


W(v) = kT, 


where k= 1.3805 X 10-23 J °K~-! is the Boltzmann constant. The factor 1/2 in (11) occurs because 
the power spectral density is defined for both positive and negative frequencies. 
3.3. The Switch 


A model for the switch discussed in section 2.1 led to eq (1). For this analysis a more general 
ys(t) composed of an arbitrary number of periodically modulated noise signals will be used. The 
voltage modulation of the ith signal [a;(t)] is analyzed into its Fourier components: 


‘e Ie 
a(t)=S Aj. med m2Tv of | 
ya 


m 


(12) 


where m takes all integer values between —© and +2. Similarly the correlation signal is 
decomposed: 


c(t)= Vi Cy,eim27vot, 
ya) 


m 
The following associated transmission coefficients are also required: 
pi(t) = a3(t), 
qi(t) = p?(t)=ai(t), 
ult) = aj(t)aj(t), 
vij(t) = a?(t)ar(t), 
which are decomposed into 


pilt ) = > P;. elmer st 


m 


qi(t)= Ss Oi, meim2mvet | 


m 


Ujjlt )= pe Ui;, medm2mvol | 


m 


vij(t )= S Vij, medmenvot : 


m 
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3.4. The y(t), Y(v), W(v) Description 


Each element of the radiometer depicted in figure 1 is described by the modification to y(t) 
before and after the element. Thus, 


ys(t)= > ailt)yi(t) (22) 
i 

ya(t)= halt) Wys(t), (23) 

yo(t)= y(t), , (24) 

yri(t) =he s(t) & yoit), 

yult)=c(t)yri(t), 

Vro(t) = Apo(t) ®& yult), 


where /ha(t), hri(t), hro(t) are the impulse response functions [20] of the amplifier, bandpass filter 
and low pass filter, respectively, and c(t) is defined in (13). 
Taking the Fourier transforms of (22) to (27) and using (12), 


Y<(v) ah Aj. — mYi(v + mv»). 


i, m 
Y,(v) =H (v)Ys(v), 
Yp(v)= Yi(v)*% Y,(v), 


Yri(v) — Avi(v)Yp(v), 


Yy(v)= > C _ mY ri1(v + mp). 


m 


Yro(v)= Art VY y(v). 


The Fourier transform of the impulse response function is known as the system function, or 
transfer function, denoted H(v), 


Hw) = | e/?7™ h(t)dt, 


and will be used to define the power response function G(v), 


G(v) = A(v)H*(v). 
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The power spectral density is obtained from (28) to (33) by means of (6) and (B-1) of appendix 
This results in 


Ws(v) — S Ai, : mé # mW i(v + MVo), (36) 


i,m ’ 
W (v)=Ga(v)W s(v), (37) 


: Li ; ; 
W v)=,., T- Yow) Yb), (38) 


Wri(v)=Gri(v) WV piv), (39) 


Wy(v)= z C_mC* ,Wri(v+mvo), (40) 


m=-2 
Wr2(v)=Gro(v)W lv), (41) 


where G,(v) (see eq 35) is the power response of the amplifier, etc. The spectral density of Wp(v) 


will not be deduced from W4(v) due to the complication arising from W4(v) being nonstationary. 


3.5. Calculation of y;2 


The calculation of yr2 (d-c output of the radiometer) is based on the following relationship 
(see appendix B): 


— __ Lim ¥Y,,(0) 
Yr2>7 >o —— 
Using (31) to (33), 
Yr2(0) = Hy2(0) S C_mH pi (mvo) Yp(mvo). 


m=—c 


The task now is to obtain Yp(myvo) in terms of the known power spectral densities of (11). The 
approach is to cast Yp(mvp) in a form to utilize the relationship (see appendix B) 


Lim Yr2(0)_ , 
roe —2—=K,, 


T 


where 


K, [’ wlwide. (45) 


and Y,,(0) is the Fourier transform of y2(t) for v=0, where y,(t) is a noise function with zero 
mean. 


If we assume a flat amplifier response so that H4(v+mvo) = H4(v), then (29) may be written, 


Ya(v) = S Ai, _pY ixa(v + pvo), (46) 


i, p 


where Yix4(v) is the Fourier transform of h4(t)#yi(t), etc. The approximation for the transfer 
function of the amplifier, H4(v) might seem unrealistic for large p, but generally Aj, , decreases 
rapidly with increasing |p|. 
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Substituting (46) into (30) will give terms of the form 
Y,(v+ pvo) WY:(v+ quo) = Yrz(v+ prot quo), (47) 


where Y,(v) is the Fourier transform of y,;(t)yz(t), x, z are any of the subscripts 14A, 2%A or 
n¥&A. Equation (47) is valid even whenx=z. For a proof of (47) see appendix B. 

We are now in a position to use (44) to facilitate the evaluation of (30). Because Y,(v+ pvo 
+qvo) gives no contribution except for x=z and v+pryo+qvo=0 (see appendix B), then by re- 
taining only nonzero terms and excluding m= 0 (because Aj, »>=0) 


’ 


Lim Yp(mvo) 7 
Too ee a é, q—mK ites; 


i,q 
where from (6), (11), and (45), for i= 1 or 2, 
: ] Pb 
Kina =5 ATi G4(v)dv. 


Thus, for the usual case that Cy = 0, (42), (43), (48) and (49) combine to give, 





Vira = NB! ‘| S CnHes(mvo) Pi, oT} (50) 


i,m 





where 


=} kHr2(0) [2G (0) G4 (0) ]*2, 


and 


i Gairndv | 


2G (0) %G4(0)’ 





or using appendix A, 
Ga(v)*%&G4(v)dv 
260) %G0)” aa 





and P;, » defined in (18) is introduced by means of (B—7) of appendix B. The term B is the con- 
volution bandwidth of the amplifier. For a rectangular power response, B is the width of the 
nonzero response in hertz (cycles per second). 


3.6. Calculation of o2 


The rms fluctuation (standard deviation) a2 is defined by 


9 
— 4 


Ths = Veolt) — Yro(t). 
The bar indicates the time average in the interval — 7/2 to 7/2. Using (10) this may be written, 
r 


x — - 
o2. = | W plv)dv — YF2- 
—s 
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If the 5-function term [23] from (38) is separated from Wr2(v) of (41), denoted by a prime, 
o3.= | W jy (v) dv. (54) 


Combining (22), (23), and (24), and assuming the modulation frequency is small compared 
with the amplifier bandwidth so the modulation factor distortion is negligible, then 


yp(t) = b ai(t)yiaale) f 
where 


Viwa(t) =ha(t) wyi(t). 


Expanding (55) and using (14) and (16), 


yp(t) =¥ pilt) ¥ixdt) + 2S) wij(t) ¥ ites (t) Y iea(t). 
i i<j 


Taking the Fourier transform of (57), then using (6) and (B—2), to obtain the power spectral density, 


, ; ' rr rrxk rr c 
W iv) = bs Pi, mP* mes, ital + MV) +4 > U ij, —mb mY ited, jxA (VT MV), (58) 


ij, 
i,m i<j, m 


where Wi, )x4 (v) is the power spectral density of yiga(t)yjx4(t) and the primes indicate 6-func- 
tion terms are omitted. It can be shown [15] that 


W 44, ices VY) = 2[(Gav)W iv) ]%& (Cay) Wiv)), (59) 
where G,,(v) is introduced in (37). Also for i ¥j, 


W snc, jxa V)= [Gav Wiv)] & (Cav) W(v)). (60) 


Using (59), (60), (10), (B—8), and (B—9) to simplify (58) and then using (39), 


W I WW) 5 ke 2 Qi, of ?+ 2 os J ij, 7'0)| G; (v)[G (vV)*G i(v)], 


i<j 


where Qj, 9 and Vij, 9 are defined via (19) and (21). 

For the computation of W y2(v) using (41), only Wyv) for v ~ 0 will contribute because of the 
narrow low pass filter F1. Consequently in using (40) and (61), Gp:(v + mvo)[Ga(v+ mvp) %&Ga(v=mr»)| 
may be replaced by Gri(mvo)[G.4(mvo)%G4(mvo)] because the power transfer functions are even in v. 
Furthermore, in the case where C, decreases rapidly with increasing |n|, then to a good approxi- 
mation G4(mvo) %®G4(mro) = G4(0)*% G0). Therefore, 


Wy (v) = ; \?M? b Qi.o T7+2 > Vij, 7:7,| Gro(v)/Gr(0), (62) 
i i<j 


where A is defined in (51) and 


1/2 
M = b & nC nGrsimvs)| 


m 
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Thus using (54), 





i 


1/2 
Or. = AMr i/2 b Qi, oT?+2 ps Vij, 7 | 
isj 





where 


sabes | ” Gralv)dv/Grs(0). 
3.7. Sensitivity 


A measure of a radiometer’s sensitivity is o71(dY¥r2/dT;). It is customary to specify instead, 
the minimum resolvable change in 7,, denoted A7T,,in. Using the criterion that a change in Yro. 
denoted Ayr2 due to a change in 7; is resolved when Ayo = oz, then 


Or2 


AT nin = aT 


Using (50) and (64), 


MS QT +25 04,00; 
i ic) 


VBr > C mH; i(mvo) Pi, n| 


m 





a a 


5 


where B, 7, C_m, Hri(v), Qi, 0,Vij,0, M are defined in (53), (65), (13), (34), (19), (21), and (63) respec- 


tively. 


a. Radiometers With Premultiplier Filter 


Most radiometers have a narrow (compared with vo) bandpass filter (denoted F1 in figure 1) 
preceding the multiplier. For this type of radiometer the sensitivity is independent of the cor- 
relation wave form. This follows from the fact that Hp;(mvo) in (50) is zero except for m=+1. 
For this case, 


1/2 
b Qi, ol? = 2 = V5, 7:0,| 
AT min oa ; Jar i 
V 2Br |P1, 1| cos W 





where 
Ci He \(vo)Pi, -1 = IC: Hpi(vo)P1, ile. 


The sensitivity for particular modulation wave forms are shown in table 1. 
b. Radiometers With No Premultiplier Filter 


To simplify the calculations, assume the transfer function of the bandpass filter is constant 
in phase and amplitude and equal to Hri(vo). In this case (67) becomes 


eis 1/2 
[c?(t)]!/? bs Qi, ol#+2 > Vy. 00| 
i<j 


AT nin = ari ere az ’ (69) 
[c(t)pi()|V Br 





where c(t)pi(t) is introduced via (B—10) of appendix B and c*(t) is introduced in an analogous manner. 
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TABLE 1. 


AT nin for special ideal switching radiometers. 


The amplifier noise is 73; Q3,0=1 for all radiometers; d is the 


fraction of the time the transmission coefficient differs from zero. 


AT min =F (Br) | S Qi.0T?+ 2S, Vis, oF; 








Modulation 


1. Square wave 
2. Square wave 


3. Voltage sinusoidal 


. Voltage sinusoidal 


. Power sinusoidal 
. Power sinusoidal 


. Power sawtooth 
. Power sawtooth 


. Rectangular 
Rect. (d < 2) 

. Rect. (d= 3) 

2. Rectangular 





Correlation 


sinusoidal 
Square wave 


sinusoidal 
square wave 


sinusoidal 
square wave 


sinusoidal 
square wave 


sinusoidal 
square wave 
square wave 
rectangular 


ae wn 


Q2, 0 Vi2,0 


Qi, 0 





nN) 


< 


1/2 
1/2 


35/128 
35/128 


3/8 
3/8 


1/3 


1/2 
1/2 
35/128 
35/128 
3/8 
3/8 
1/3 
1/3 
l—d 
l—d 
l—d 


YN NY 
< 
wo) 


S 


yw 





al(V2 sin 7rd) 
1 


(l—d)>" 























[d(1—d)]- +2 


lied 


*Results of author conflict with this paper as noted in footnote. 

°Results of author agrees in detail with this paper only when 7, = 72 

+The calculation made was equivalent to assuming 7; =0 °K. 

+A rectangular filter bandpass assumed in calculation. 

® Dicke [3] obtains a sensitivity greater by a factor 2 (thus AT; less by factor 2) than this paper. 
» Goldstein [4, 5], a minor algebraic error corrected in ref. j. 

© Bunkin, Karlov [17]. _ 

4 Strom [6, 7] obtains a sensitivity less by a factor V2 than this paper 

e Strum [8]. 


‘ Kelly, Lyons, Root [11, 12] agree only when radiometer balanced. 


The bandpass defined does not agree with this paper except for rectangular bandpass. 


* Colvin [13] agrees only for balanced, constant power radiometers. 
” Knight [9]. 


occur. 


Equations (3) and (4) of Knight's paper omit a factor 27. In addition, Knight's results differ from this paper in that the terms in 7; and T2 do 


Ward and Richardson [18] obtained a sensitivity greater by a factor 2 V2 than this paper 


’ Johnson [10] algebraic error for square wave modulation. 


* Tiuri (14). 


4. Special Results 
A variety of particular radiometers have been considered in the literature (see sec. 2). These 
radiometers were listed in table 1. The radiometer types shown in table 1] are ideal in three ways. 
First, the switching element introduces no noise of its own, second, the maximum transmission 
coefficient is unity, and last, the relative phases of switching and correlation wave forms a(t), b(t), 
and c(t) are selected to minimize ATijn. The factor F is equal to 2-"/?|P,,,|~! for sinusoidal cor- 
relation or for radiometers with a narrow band precorrelator filter (and any reasonable correlation). 
These radiometers are denoted as sinusoidal correlation in table 1. For other radiometers, F is 


equal to [c%(t)]"[e(t)pi(t)}-!. 





Explicitly the power modulation wave forms used are 
(i) square wave modulation, 


| ieee 
- ¢ 


Ss cumanagmpaians eim2m ef 
mir 
x 


pi(t) = 
m=— 
(ii) voltage sinusoidal modulation, 


pi(t) =3/8+ 1/2 cos 27 vet + 1/8 cos 42 vot; 
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(iii) power sinusoidal modulation, 
pi(t) =1/2+1/2 cos 27 vot: 
(iv) power sawtooth modulation, 


Pilt) =12+j S (2ark)-reik?nv0!; 


k#0 


(v) rectangular modulation, 


sin kid . 
mls) = > rw’; (74) 
2 ki 


where d is the fraction of the time that p,(t) differs from zero. Except for voltage sinusoidal 
modulation, p2(t) is chosen such that p:(t) +p2(t) =1. 
4.1. Nonideal Switch 


The theory developed applies to arbitrary switching wave forms. There are two common 
deviations from the ideal forms used for table 1. First, the nonideal power transmission coefficient, 
pi(t) of the switch varies from zero to Pi, max instead of from zero to one. Thus. 


pilt) =i, maxDilt) 


so that 


Q; o= Pt. my 05 


and 


Pe hs ats ; 7 wi 
J ij,0 Pi, maxPj, max ij, 0- (76) 


Equations (75) and (76) may be used in place of the corresponding ideal parameters, to extend the 
results of table 1 to the switch with nonideal transmission coefficient. 

A second complication is that the nonideal switch radiates thermal energy in direct proportion 
to its temperature, and the fraction of incident energy absorbed. For a matched three-port device 
at temperature 7,,,, the effective temperature of the radiation p.y(t)T.y is [24] 


Psw(t) T sw = [1 — pilt) — pelt) |T sw, (77) 


where the p(t) are given by (14). The thermal noise from the switch may be treated formally as a 
noise source of the temperature of the switch which is power modulated by amplitude p,,(t). Thus 
for the ideal square-wave modulation, rectangular modulation, power-sinusoidal modulation, and 
power-sawtooth modulation of table 1, 


Pew (t) =0. (78) 


For ideal voltage sinusoidal modulation, 


Psw(t) = cos 47rvot. 


1 
4 
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TABLE 2. Balanced, constant switched power, nonideal radiometer 
sensitivities. 
The maximum power transmission coefficient through the switch is 
Pi. max= P2, max=@: T:=T2: the effective amplifier noise is Th: 
and the physical temperature of the switch is T.y. 


Ti: +a7'T,+ (a7! —-1)T sy 
|? 
VBr 





AT nin ae F 


A. Narrow band precorrelator filter or sinusoidal correlation 


Modulation 


Square wave 
Power sinusoidal 
» 5 ‘ 

Power sawtooth 


Rectangular* (V2 sin 7d) 
B. No precorrelator filter 
Modulation Correlation 


Square wave square wave 
Power sinusoidal square wave 1 

Power sawtooth square wave 4 
Rectangular* rectangular [d(1—d)]-!? 
Rectangular* (d < >) | square a” 
Rectangular* (d > 3) | square (i—@* 








*d is the fraction of the time the transmission coefhcient is dif 
ferent than zero. 


4.2. Balanced Radiometers With Nonideal Switches 


The expression for A7,\, simplifies when the noise sources being compared have nearly the 
same effective temperature and when the power to the amplifier is constant [i.e., p(t) +po(t) 
=constant]. Table 2 lists the results for radiometers having nonideal switches that satisfy these 
special conditions. The sensitivity indicated is for the correlation phase which minimizes AT pip. 

The AT), for the voltage sinusoidally modulated radiometer does not simplify for 7; =Te. 
However, for a radiometer with a narrow band precorrelator filter (or sinusoidal correlation) and 
a switch where the maximum power transmission coefficients pi, max = P2, max = @: AT min for T:=T2 


is bounded by the following expressions: 


3 7 
E Ti, +a Ts +(a oe Tse | 
AT min > 2 v2 > 
Br 





‘ 


g lita 'T;+(a 1—— 





AT nin < 2 v2 | 


~ . ° [> - 
For no precorrelator filter and square wave correlation, 2 V2 is replaced by 7. 
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5. Experimental Procedure 


The sensitivity of a Dicke radiometer similar to Wells, Daywitt, and Miller [1] was deter- 
mined experimentally. A switchable circulator was used for the switch instead of a magic tee 
and motor driven attenuators. The switchable circulator provided square-wave modulation at 
1000 Hz. The 30 MHz output of the IF amplifier was attenuated by various amounts and then fed 
to a crystal detector. The detector law of this detector is shown in figure 2. The arrows indicate 
operating points for which the radiometer sensitivity was measured. The numbers near the 
arrows are the minimum resolvable temperature changes measured for each level. 

The mixer-preamplifier-IF amplifier voltage response is shown for one channel in figure 3. 
The other channel is similar. The markers along the bottom of the figure are synthesized at 2 
MHz intervals. The convolution bandwidth was calculated using a computer program based on 
(53). The system bandwidth is determined from a single channel response as discussed in ap- 
pendix C. 

The time constant for the low pass filter of the lock-in amplifier nominally set at 10 sec inte- 
grating time is based on the voltagetresponse measurement shown in figure 4. ‘ 

The recorded data were plotted on probability paper. The ordinate of this paper is the am- 
plitude (voltage) of the radiometer output, and the abscissa is the percent of the time the recorder 
trace exceeds the ordinate value. The scale of the paper makes the plot a straight line if the 
data have normal (Gaussian) distribution. From the plot, the mean of the data is the value of the 
ordinate for which the abscissa is 50 percent. The standard deviation o (for normal distribution) 
is the difference in the ordinate values for abscissa values of 50 percent and 84.1 percent (or 15.9% 
and 50%). Figure 5 shows two sets of data, one before and one after a change of 0.0087 dB (21.6 
°K change from 10,800 °K) of input power is made in one noise source. 


The various parameters for the radiometer under test are given in table 3. The theoretical 
ATmin is based on table 2 for the balanced radiometer using square-wave modulation. The 


for a AT, = 21.6 °K. 

A feature common to all data and shown in figure 5 is an increase in o when the noise power 
decreases from null. No change in o is noticeable when the noise power increases (out to 0.2 
dB change). The cause of this feature is not known. 

The radiometer sensitivity as noted in figure 2 for the radiometer at various operating points 
are comparable. 


6. Conclusions 


The areas of agreement within the literature have been ferreted out and the nature of dis- 
crepancies determined. A new analysis of the Dicke radiometer was presented which was general 
enough to include all of the papers reviewed as a special case and to extend the sensitivity calcu- 
lations to include radiometers with a lossy switch. 

The most sensitive radiometer is a square-wave modulated radiometer with wide band (i.e., 
no premultiplier filter) square-wave correlator. However, the potential 11 percent advantage in 
sensitivity for wide band square-wave correlation over systems using a narrow premultiplier filter 
must be weighed against the practical difficulties associated with wide band systems. Also the 
narrow band correlator insures that the output of the radiometer is proportional to the difference 
of input powers [i.e., to T;—T2; see eq (50)] independent of the asymmetry of the modulation. 

The sensitivity of a square-wave modulated radiometer using a narrow precorrelator filter was 
determined experimentally. The measured sensitivity was the same as the theoretical radiometer 
using a full-wave, square-law detector. The sensitivity of the radiometer was independent (within 
the experimental error) of whether a half-wave square-law, linear, or intermediate-law detector 
was used. 
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FIGURE 3. Amplifier voltage response. 


TABLE 3. Parameters of the Dicke radiometer under test 


Predetection bandwidth (B) 9.2 MHz+ 10% 
Postdetection integrating time (7) 18 sec + 10% 
Receiver noise temperature (T7',) 750 °K +5% 
Maximum switch transmission (a) (2.57)-' 4% 


Switch temperature (TJ sw) 295 °K 

Noise sources (7; = T2) 10,800 + 2.5% 
Switching frequency (vo) 1000 Hz 
Switching waveform Square wave 
Microwave frequency (local oscillator) 9 GHz 

AT min (theory) 2.3 °K + 15% 
AT nin (experimental) 2.3 °“K + 20% 
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7. Appendix A. Equivalent Statements for B and 7 and Notation Conversion 
for Literature 


Various statements are used in the literature for B and r. This section gives the definition 
used in this paper in alternate forms which are easier to compare with the literature. 

Using the Wiener-Khintchine Theorem, (7), and assuming the order of integration can be 
interchanged, the following relationships can be proven: 


G(0) #%G(0) =| GXv)dv, 
[’ Gv) *® G(v)dv = ih G( vide 
g(0) =| G(v)dv, 


Ho) =|" h(t)dt, 


and Parsevals relation, 


[’ Hw Hid = | h?(t)dt, (A-5) 


where g(t) is the Fourier transform (autocorrelation function) of the power response G(v), and 
h(t) and A(v) are a Fourier transform pair, with h(t) a real function. 
If instead of using frequency v in hertz (cycles per second) as a parameter, the angular fre- 
quency w= 27v is used, then the following conversions may be used: 
v=a/(277), 

Gv) = 27G.(o), 

W (v)=27rW.(o), 

Y,(v) = 27rY..(@), 

H,(v)=27H.\o), 

GAY) ®& Gv) = 27G,(@)*&G.(o), 


where the subscripts v and w are introduced to distinguish a difference in functional form. This 
subscript will now be dropped. The functional form appropriate to the function used is under- 
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stood. The conversion from W(v) to W(w) is chosen so that 


W(v)dv = W(w)dw. 


In the » notation, 


e Joly t)dt, 


Lim 


W (w)= 1 (27/T)Yiw)Y*(w). 


The conversion system above is not unique. The easiest way to determine what system a given 
author is using is by examining the form used for the Wiener-Khintchine theorem, or statements 
which include the autocorrelation function. The autocorrelation function is universally defined 
in a manner equivalent to (9). 


7.1. Equivalent Statements for B 


ih Garde] 


26.4(0) %&G4(0) 


The definition of B is in (52): 





Using (A-2), 


[’ Galv)*%G4(v)dv 
aia 26.(0)%G4(0)  ” 





or (A-1), 





| ; Gade | 
B=—= 
2 | G?(v)dv 


Beginning with (A-8), noting that G4(v) = H4(v)H}f(v), and converting to @ notation, 


oi HHao)/*do] 


an | |H4(w)|4*dw 





Beginning with (A-8), applying (A-3) to the numerator and (A-5) to the denominator, 


B= g*(0) 


or a (A-10) 
2| gi (t)dt 
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Using a normalized autocorrelation function, pa(t)= ga(t)/g4(0), 
2B)-1=| ~ pX(t)dt. 


7.2. Equivalent Statements for 7 


The definition of 7 is from (65): 


. = {7 Gro(v)dv/Gr2(0). 


Noting that Gro(v) = Hro(v)H*.(v) and converting to the w convention, 
to) F2 2S 


[’ |H-2(w)|2dw 
1 = ° 





277 H?..(@ = 0) 


Using (A-4), and (A-5), and assuming a causal system such that A(t) =0 for t < 0, 


| h2..(t)dt 


ih hystthde| 


7.3. Literature Conversion 
Several representative statements from the literature referring to the minimum resolvable 
change in temperature for a radiometer near balance (7; = 72) and for AT pin < Tn will be converted 
to the notation of this paper. 
a. Dicke [1946] 
Dicke’s [3] equation (21) for square-wave modulation and cosine correlation is, 


9 


[Jiri] sa 


sio) |” F (w)|?dw 





Using (A-9) and (A-13) and the following conversions, 


AT= AT nin. 


N =(T,+17;)/7; (for 7; =290 °K), 
F(w)=H4(o), 


S(w)= H; 2(@), 


and (A-15) becomes 


> 


AT nin = a(T; + T»)(8Br)- 2. 
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b. Goldstein [1955, 1956] 


For voltage sinusoidal modulation and sinusoidal correlation Goldstein [4, 5] obtains 
40°, V y/a. 
Goldstein’s work corresponds to 7;=T72=0 °K and 


Of, = KAT mins 


where / is Boltzman’s constant. Thus (A—16) becomes 


AT min = 40 n(2Br)- 1/2. 


c. Bunkin and Karlov {1955} 


For square-wave modulation and sinusoidal correlation, Bunkin and Karlov obtain [17] 





T T2 1/2 
AAD AT (1 AT }}" oe 


T= a Tn TG 30, a0 |! TF 3 Tet To 


If the precorrelator bandwidth AQ, is much greater than the postcorrelator low pass filter band- « 


width AQ, then 


AQ, AQ 


AQ, + AQ =a. 


Using the following conversions 
87 = AT nin. 
Tin me Tn, 


Te=T;, 


for AT < 1, (A-17) becomes 
AT... —_ aT; - T,)(2Br) 1 2. 


d. Kelly, Lyons, and Root [1958, 1963] 


For square-wave modulation and detection, (49) and (31) of Kelly, Lyons, and Root [12] becomes 


+T (0) | (A-18) 


AT=2T,+T,) fo 
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Equations (A—11) and (A-14) are used to establish the conversion for B and rt. 


Using the following 
conversions, 


(A-18) becomes 


AF wan = 2 T; ar; i i \(Br)- 2, 


e. Knight [1962] 


For rectangular modulation and rectangular correlation, Knight [9] obtains 


, ' all , 
S_ 27 a(o}%—o%) a (1 a ) 
N 21loa% a+b at+b/’ 





where the 27 omitted in Knight’s eq (4) is added. 


The appropriate resolution criterion is S/V=1. 
With this condition, then 


= <r 
o2—02=hAT mins 


and (A-19) becomes 


AT nin = Tn{ Brd(1 — d)]- *2 


8. Appendix B. Theorems on the Fourier Transtorm 


Several theorems involving the Fourier transform of the time function y,(t) and y-(t) of zero 
mean are proved in this appendix: 


Y(v)¥*(v) _ 


Lim : 0, 
Tx 7 





for yz(t) uncorrelated with y(t), and 


2 Cee 
Lim 7 Yec(0)¥*2e(v') =0, 


unless yz(t) = y(t), yz(t)= yi (t), and v=v’, and 


Y(0) 
T° 


y= Lin 


T-« 


Lim a a W (v)dv, 


Y,(v + pvo) ®& YAv + qvo)= Yr(v + pvo + quo). 


Pi,, 


] (° 2 
a3(t)e?™"oldt, 


~ OJ e/2 


~ * 
Qi.o = > Pint tos 
m 


Vio= dl ij, m Uj (B-9) 


ij,m-? 
m 


where P;,, and Aj, p are defined in (18) and (12) respectively, and 0 = v7 ' is one period of the periodic 
signal p(t), and 


SY CaP, a=ct)pilt) (B-10) 


where 


niscaiaiaiaaatancii ] 8/2 
cpl = 5 | c(t)pi(t)dt. (B-11) 


0 


0/2 


8.1. Theorem 1 


Taking the inverse transform (3) of the left-hand side of (B-1), then taking inverse transform of 
Y.i(v) and Y2(v) and interchanging the order of integration, one obtains (Lim T — & is understood) 


T | e2™Y(v)YF(v)\dv=T aT yu(')yat’)at’ae” efter + Ody, 
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The integration over v yields the Dirac delta function 5(t—t' + t”) [24]. 


Thus integrating over t’, 


T tf eam WW o)dy=T f° yilt +t"yo(t")dt”. 


If yi(t) and y2(t) are not correlated, then the integral on the right is finite, and when divided by T 
tends to zero as T>%. Thus, the left-hand side is zero, the Fourier transform is also zero and 
(B-1) is proved. 

8.2. Theorem 2 


The theorem is proved in a manner analogous to theorem 1. 


8.3. Theorem 3 


By definition, 


= . 1 (tl 
y=Lim= y(t)dt. 
To r/2 


|>T/2. Thus, 


The function y(t) is truncated so that y(t) equals zero for |t 


y(t)dt. 


x 


= : ] 
y=Lim= 
| edi I 4 
Taking the Fourier transform of y(t) and interchanging the order of integration, 


; ; ] 00 % 
y= Lim T Yio) | er de. 


To 


The integration over time yields the 5-function of v. Subsequent integration over v yields (B-3). 


8.4. Theorem 4 


Let Y,,(v) be the Fourier transform of y,(t)y,(t) where y,(t) and y,(t) are random variables with 


zero mean: 


x 
d = Qnty, ’ 
Yx2(v) | e~ 2avty (thy, (t)dt. 
If y,(t) is not correlated to y,(t) or e~'?™, then the mean of the product is zero. Thus, Yx,(v) is 


finite and if divided by T, will approach zero as T—>~. If y,(t) = yx(t), then 


Te 
Yxtv)= | e~ 27ty2(t)dt, 
-T/2 


where the limits of the integral are taken from —T/2 to T/2._ This does not change the value of the 
integral because y(t)=0 for |t}=>T/2. When y2(t) is rewritten as y2 +[y2(t)— y2] and the above 


equation is divided by T (the limit T > © being understood), 


ey i 6: T/2 sated 
T~1Y,,()=T ye e~ mtd +| e ara ty2tt)— veld} 
_T/2 _ 7/2 
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The integral in the second term is finite and the term goes to zero as T—>. The first term on the 
right-hand side goes to zero unless v=0. For v=0 the term approaches yz. Thus, 


Saal y2 for v=0 and x=z, 
Lim Y,,(v)=47 
Tox 0 otherwise. 


Using (10) yields (B-4). 


8.5. Theorem 5 


Using the definition of convolution, eq (5), 
x 
y 7 4 , , ; 
Y,(v+ myo) ®Y (v + qvo) = Y,(v’ + mv) Y,(v — v’ + qvo)dv’. 
x 
Taking the Fourier transform and interchanging the order of integrations, 


avol'ly (t)yz(t')dt dt’ - eo ea" dy’, 


J x 


x 
y x(Vt mvVo)® Y avt qvo)= [| e j27[Mvot + vt’ + 
x 


The integration over v' yields &(t—t'), then the result of integrating over t’ is identified as 
Yx2(v + mvp + qvo). 


8.6. Theorem 6 


To obtain (B-6), both sides of (18) are multiplied by e*27"’*' and then integrated over a period, 
using for p;(t) the equivalent a(t) as noted in (14). 


8.7. Theorem 7 
Using (B-6) and (12), 


l me j on myierd 
Pir=6G >| , A genrnrr FA, emaeret id. 


m,s 


Due to the orthogonality of different Fourier components, the integral has nonzero value only when 
s=r—m. This leads directly to (B-7). 


8.8. Theorems 8, 9, 10, and 11 


The proof of (B-8), (B-9), and (B-10) are analogous to the proof of (B-7). 


Equations (13), 
through (21) are needed for various definitions. 


9. Appendix C. Radiometer Using an IF Amplifier 


A schematic for the IF amplifier is shown in figure C-1. 


Instead of the amplifier shown in 
figure 1, now a mixer, local oscillator, and IF amplifier are used. 


The amplifier noise temperature 
T,, of figure 1 is now interpreted as the noise temperature of the mixer-IF amplifier as referred to 
the input of the mixer. Using the same notation convention as the body of the paper, 


Yio(t)=7 cos 274 t, 
y¥sm(t)=nys(t) cos 27vit, 


ya(t)=ha(t)¥ysmu(t), 
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where 7 is the amplitude of the local oscillator of frequency vi. From (C-1) 


, it is clear that an 
ideal balanced mixer is assumed. Taking the Fourier transform, 


Ysu(v)= 2 [¥s(v t+ v1) + ¥s(v—v1)], (C-2) 
Y s(v) =H a(v)Y¥ sy(v). (C-3) 


Applying (6), 


W sulv) = (77/4 s(v + v1) + Ws(v—v1)], 


W 4(v)=C4a(v)W su(v). 


9.1. Calculation of V2 


The calculation of ¥rz proceeds in a manner similar to section 3.5 except (C—2) and (C-3) are 
used in place of (29). Thus, instead of (46), now we have, 


Y,(v)= 2 > Ai, plYixa(v + pvo + 1) + Ying a(v + pvo— 1) ]. 
i,p 


Substituting (C—6) into (30) will give terms of the type [see eq (47)], 


Y-(v + mvp + 11) & Y2(v + quo + 11) = Yre(v + mvo + qvo) 


and 


Y2(v + mvp + 11) ®& Y(v + quot 1) = Yrz(v + mvo + quo + 21). (C-7) 
Using an argument similar to that used to obtain (48), noting that the terms like (C-7) do not 
contribute, 
‘ moti ] , . 
Lim ee ~ 97 es i, q4 ~q mKita- (C-8) 


Equation (C-8) is identical to (48) except for the multiplicative constant } 7. 
of the c 


Thus the remainder 
iyi mee wie ' nee Tees: 
‘alculation is unchanged except that the gain functions A of (51) now contain 3 n?. Thus 


= 5 wkHr2(0)[264(0) & E4(0)]"”. (C-9) 
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9.2. Calculation of o> 


The calculation of or: proceeds in a manner similar to section 3.5 except (C—4) and (C-—5) are 
used in place of (37). Because Ws(v + v1)=Ws(v — 1), (C—-4) and (C—-5) combined become, 


W 4(v) = 3 WGa(V)Ws(v). 


Thus the calculation is similar except Ws(v) is replaced by } y?Ws(v). If (C-9) is used for the 
definition of A, then the statement of yr2 and oz in (50) and (64) are otherwise unchanged. 


9.3. Comments 
The bandwidth B is defined in a similar manner whether the amplifier is at the signal frequency 
or an IF frequency. For an IF amplifying system, B is the bandwidth of the IF amplifier. If for 
some reason the signal were only in the signal channel of the IF amplifier (the noise remains in 
both channels), then 7; and 72 would be replaced by 27; and 272, but no other changes. Thus, 
the single channel IF radiometer would be equivalent (as far as sensitivity is concerned) to a similar 
double-channel IF radiometer except it would have twice the amplifier noise temperature. 
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lhe correlation radiometer is analyzed to determine the sensitivity that can be obtained under 
various operating conditions. 

The radiometer using a sine wave comparison signal is analyzed and compared with the usual 
radiometer that employs a random noise for the comparison signal. I[t is found that the radiometer 
employing the sine wave comparison signal is the more sensitive of the two circuits, particularly in 
the case that the effective temperature of the input noise signal is greater than the effective input 
temperature of the amplifiers. 

It is shown that if nonidentical amplifiers are used in the correlation circuit, the properties of the 
radiometer are determined by the portion of the amplifier response functions in the frequency interval 
that the two response functions overlap. The effect of amplifier gain fluctuations are considered, and 
although the correlation scheme reduces the effect of gain fluctuations, it is shown that they still do 
contribute to the output fluctuations of the radiometer. 

Calculations are included showing that the effect of a differential phase shift between the two chan- 
nels is a reduction in radiometer sensitivity. The same conclusion is reached concerning the effect of a 
differential time delay. 

Finally. it is shown that if the comparison signal and the input signal have the same statistical 
properties, the requirements on the multiplier are less stringent than if the two signals have different 
statistical properties 


Key Words: Correlation radiometer, differential time delay, gain fluctuations, imperfect mullti- 
plier, noise comparison, nonidentical amplifiers, sine-wave comparison signal. 


1. Introduction 


The use of correlation techniques in radiometry has been suggested recently by several workers 
[Strum, 1958: Blum, 1959: Colvin, 1961; Allred, 1962]. One of the major reasons for this interest 
is that by using correiation techniques, it is possible to build a radiometer in which amplifier gain 
fluctuations contribute less to the fluctuations present in the output of the instrument than they do 
in the conventional radiometer: thus, the correlation radiometer might be expected to be superior 
to the conventional radiometer under conditions in which gain fluctuations are important. Strum 
[1958] has pointed out that as the noise level in a conventional radiometer system is reduced, the 
fluctuations due to changes in amplifier gain become more serious: therefore, at cryogenic tempera- 
tures it is very important to reduce this contribution to the total radiometer fluctuation. Thus, 
it was felt to be desirable to calculate the sensitivity obtainable from instruments of this type. 

The analysis was carried out on the correlation receiver shown in figure 1. Although several 
different types of correlation receivers have been suggested in the literature, the operation of all 
of them is quite similar to the one chosen here. Therefore, the results obtained here should be 
applicable with appropriate modifications to many of the correlation radiometer circuits in the 
literature. 

The.first problem to be considered (see sec. 2) is the calculation of the sensitivity of a radiom- 
eter with amplifiers whose gain and passband are identical but with different effective tempera- 

*This research was supported by the Advanced Research Projects Agency under Project DEFENDER and was monitored by M. |. Witow under ARPA Order 
No. 515 


**Radio Standards Laboratory, NBS Laboratories, Boulder. Colo. 80302. Permanent address, Colo. State University, Ft. Collins, olo. 80521 
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tures at the input circuit. These calculations are carried out for two types of input signals: 
(1) both input signals consist of white noise and, (2) one input signal is white noise and the other 
is sinusoidal. The first case was selected because it corresponds to the problem considered by 
Colvin [1961], and the second was chosen because it was the problem considered by Allred [1962]. 
Throughout this calculation, all filters are assumed to have a square passband. This assumption 
is made because it allows the results of the calculation to be written rather explicitly and thus it 
is easy to see the physical significance of the results. When it is desired to apply these results 
to a real radiometer with filters that are not square, it is only necessary to use the formalism of 
Colvin [1959] that expresses the effective bandwidth of a filter in terms of a convolution integral. 
Thus, this generalization can be carried out rather simply. 

Section 3 is an analysis of the correlation radiometer in which it is assumed that the two ampli- 
fiers have different gains and different response curves. 

Section 4 deals with the same circuit, but here the effect of variation in phase and in time delay 
of the signals in the two amplifiers is considered. 

Finally, section 5 deals with some effects produced by an imperfect multiplier. 


2. Sensitivity of an Ideal Correlation Radiometer 


The circuit to be analyzed is shown in figure 1. Two noise sources, designated X and Y, are 
connected to opposite arms of a matched hybrid junction. The signals from the other pair of arms 
are applied to amplifiers whose transfer functions are R,(w) and R2(w) respectively. The outputs 
of the amplifiers are multiplied together and the product is filtered with a smoothing filter (low 
pass filter). The output of the smoothing filter is displayed on a d-c instrument. 

As the following analysis will show, the average value of the deflection of the output instru- 
ment is proportional to the difference in noise temperature of the two sources, X and Y. Thus, 
when the two noise sources have the same effective temperature, the average value of the output 
deflection is zero and the radiometer is said to be balanced. The instantaneous deflection is a 
stochastic function of the time; thus the random fluctuations of the output deflection produce an 
uncertainty in the experimental conditions that correspond to a balance. 

The sensitivity of a radiometer is usually defined as the change in temperature of one of the 
noise sources that will produce a deflection whose magnitude is equal to the root mean square of 
the output fluctuation. Therefore, in order to calculate the sensitivity of a radiometer it is neces- 
sary to obtain expressions for both the rms value of the output fluctuation and the change in aver- 
age deflection per degree change in effective temperature of one of the noise sources. 

In this section the assumption will be made that the two amplifiers have the same power trans- 
fer function, and for simplicity, a “‘square passband” will be assumed. ' 

If x(t) is the signal voltage due to source X that leaves each amplifier, and y(t) is the signal 
voltage due to source Y that leaves amplifier 1, the signal voltage — y(t) due to source Y must leave 
amplifier 2. It will be further assumed that the amplifiers introduce noise voltages z; and z2 
respectively. 

Since the random functions x(t), y(t), z:(t), and zo(t) represent noise from physically separate 
noise generators, it will be assumed that they are mutually independent. It is also assumed that 
each is a second order, stationary Gaussian process possessing a continuous spectral density 
function. These assumptions imply that the processes are ergodic so that:time averages and 
ensemble averages coincide. It is also assumed that each of these functions has zero mean. 
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If the output of amplifier 1 is designated by v; and the output of amplifier 2 by vz, then 


we—xtytz 


The multiplier will be assumed to have a response law such that 
U= avjd2 (2) 


where u is the voltage output and a is a constant of proportionality. 

For any random quantity é, let € be the expectation or ensemble average of €. (This is the 
equivalent of the notation E(é) used in the literature of mathematical statistics.) Using the inde- 
pendence and normality of the random processes and the fact that each process has zero mean, 
the autocorrelation function of wu is 








WT) = u(t u(t + 7) =a? vi(t)volt v(t + r)v2(t +7) 





= a? [x(t)— y(t) + z(t) (x(t) — y(t) + 22(t) (x(t) + Mt) + z(t )zolt)| [x(t + 7) — (t+7)+. 
WT) = a?[ rl) + Wy) — WW AD KO) + We (7) (WAT) + WAT) + We(7)(WA7) + Wf 7)) + We (T)wz{7)]. (3) 


If €;, i=1, 2, 3, 4 are Gaussian variables, each with zero mean and 
figj = Oi 


then 


E1€2€3&4 = 012034 + 013024 + 014023. 
Setting €;=&=2(t) and €;= &,=2(t+7), there results 


Walt) = x7(t)x%(t + 7)=22(7) + W2(0). 


When (4) is applied to (3), the result is 
a~*Wd 7) = 2[UR(7) + YR(z)] + [{0) — ¥,(0)]? + We,(7) [W(7) + wf7) 
+ we) [Wr(7) + bt) ] + We, (T)b2{7). (6) 
The Fourier transform of the autocorrelation function of a process is the spectral density 
of the process. Further, the Fourier transform of a product of functions is the convolution of the 
Fourier transforms. Thus, designating the spectral density by Q, and the convolution operation 
by *, | 


a~*Q.A f) [ Ws(0) = W,(0) |25(/) +2 [O, * VAf) =a Qy * Q0,(f)] 


+ [Qz,+ Qe] * [Q2+ Qy](f) + Qz, * Qz(/). (7) 
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Case 1. Consider both X and Y to be white noise generators at temperature 7; and Ty respec- 
tively. When the spectral densities are evaluated, it is necessary to consider the power division 
that takes place at the hybrid junction. These spectral densities are (see appendix) 


0.40, ; (kT ,)PB, (8) 


Ox *Qz, ~ ; kT ,T.,B, (8a) 


O:,* Oe ~ 5 KT: T2,B. (8b) 


Thus, using Q’ to designate the portion of the spectral density that contributes to fluctuations, 
this portion of (7) becomes 


a*k?B 


Qu=F— (T3413 + THT, + Te.) + Ty(Ts, + Taq) + 27 2,7 ey]. (9) 


With the assumption that the two amplifiers have the same effective noise temperature, this 
becomes 


ah?B 


= 
aaa 


[T3+ T7+ 27,T,+ 2T,T, + 272). (10) 


If the smoothing filter has a power response G(/), the output power spectrum from the smoothing 
filter is 


W( fp =GNQu/). (11) 


It is well known that integrating a power spectrum over all frequencies results in the mean square 
of the output voltage or current; so if the output voltage from the smoothing filter is w(t), and 
G(f) represents a square pass filter of width 6, then 


a a (df= : a®k*GobB [T2. + T2 + 2T2 + (TT: + T,T>)] (12) 


where Go is the power response of the filter at zero frequency. The fluctuations at balance are 
required, so the balance condition is substituted into w?. The balance occurs when 7,=T,, so 
for balance 
uw? = ak?’GobB(T + Tz). (13) 
Therefore, in terms of an rms voltage 
Wems = (w2)"? = ak(GobB)(T, + 72). (14) 
Next, we calculate the deflection arising from a certain imbalance of the input signals. From (2) 
U = aviV2 = A(x + Y + 21) (X — Y FT Za), 
= alx® — y*? + 2122 + 21(x — y) + 20o(x + y)]. (15) 


Since x(t), yt), z:(¢) and z(t) are all uncorrelated, when (15) is averaged over time, the result is 


u =a(x?— yy). : (16) 
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Again applying the theorem that the integral of the power spectrum over all frequencies is the 
average square of a function, (16) becomes 


u=a| | oud—| oar 


kB 7 
aoe eT. 


In order to examine the effect of a small imbalance, set 
T,=T,+ AT. 
Then 


u+Au= ant [T—(T-+ AT). 


Since u=0 when T,=T,, 


|Au| =“22 ar. 


The resulting signal at the output of the smoothing filter is 


|Aw|= Gy? |Aul=< GaekBAT. (20) 


As is usual in radiometer calculations, the assumption is made that the minimum detectable signal 
occurs when the deflection is equal to w,m,. Thus, equating the right-hand sides of (14) and (20), 
the condition corresponding to minimum detectable signal results; i.e., 


ak(GobB)*/>( fiz + T:) = 


< GV2kBAT. (21) 


Thus, the minimum detectable temperature difference is 
AT = 2(b/B)'(T, + T:). 
Case Il, Let Y produce a sinusoidal signal, so that 


QOAf)= 4 


The expression for Wy2-y2 is needed, so it will be calculated first. From the previous calcu- 
lation, it is known that 


Ware a Ware + Wye io 2x” . 


Thus, it is necessary to compute W,p: 


Wye = V7(t)y(t +7), 





,4 
u é ‘ 
=— cos? wot COs? Wo(t +7), 


4 


a = 
— cos? wot +2 sin? @oT |: 


8 
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Making use of the above and (5), the result is 


ieee sae 
i. 


3 ] 
Wre- pax + Ay.) + 2 m ah cos? WoT + 5 sin? wor) 2x “@ 


From this, the power spectral density is 


On-n=| Wr2—ye "dr, 


=(# -#2) Bf) + 2(Qz* Q)+7 ah a(/) +22 af -2fs) += (+ 2fa)| 
Therefore 


Qu=a@?Qr,r. = “| Qe yz + Qz(2-y) + Qelrty) + On 


me -*) an+z 7 [af f;) + Bf + 2fo)] + 20x * 


+ Qz, *(Qzr + Qy) + Qz, *(Qzr + Qy) + Qz, * Veet" (29) 


Since the present calculation is to compute the fluctuations occurring in this radiometer, only 
the terms that contribute to the fluctuation need be considered; therefore they will be expressed as 


Qi = @{2(Qr*Qz) + (Qz, + Qz.)*(Qz + Oy) + Oz,*Qz.}, (30) 
where Q.(/) denotes the fluctuating portion of Quy. 


All of the evaluations of the convolutions that occur in (30) can be taken from (8), except 


Q,*Q:,; which is, 


Q*0-,=0:.40y= | Q2f QWf—-f df’. 


where Q, is given in (23). Thus, 


Qvt0e,= | Qalf) F (BUS +4) + US — fd 


v* (31) 
=7 10:,°+fo)+ 02, fo). 


Since Q., is assumed to be constant through the passband, when fp is in the passband and f is 
small, this convolution becomes 


Qy*Qz, = 


Therefore 


’ 9 ] 9 ] 9 6 y ‘ 
O,=@ iq (kT 2)? +4 PBT AT, + T2,) +1 py. +T.) +4 eT.,T.B}- 


4 2 
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By analogy with ihe previous calculation, the spectral density of the output of the filter can be 
written as 


W(f)=GCNQAf), 


and 


w= | ” Wfdf=a2Go(2b) F (KT 2B +4 BT ATs, + Tx) +4 vP(Ts, + Tx) + 5 


eT.,T \. (34) 


The d-c term of (29) indicates that the system is balanced when k7,B=v?/2. Also, to simplify 
the expression, it will be assumed that T.,=T:,. Then, at balance 


w == aGobk?B (5 7 E+ 27.7 v T:) . (35) 


By analogy with the previous calculation, a root mean square fluctuation amplitude may be de- 
fined by 


1/2 
Wma = a GobB)"2 (5 T3+20.T:+T2) (36) 


As mentioned previously, in order to arrive at the radiometer sensitivity, two parameters 
must be computed. The first is w,ms, which is given in (36). The second is the change in deflec- 
tion due to a small change in temperature of the thermal source. The deflection sensitivity will 
be the same as it was in the previous calculation, so it is 


‘“ akB 


|Aw|= Gi? —— AT: 


Again, the minimum change in temperature that is observable will be taken to be a change 
that produces a deflection equal to the wyms. Therefore 


AT = alKGoB)" (3 T2+27,T 


ae akB 


0 9 


and the final result is 


ee. ee = 
aT=2 (5) (5 T$+27.7.+T) (37) 


When this result is compared with (22), the sensitivity of a radiometer with both signals as- 
sumed to be white noise, it is evident that this instrument has a somewhat higher sensitivity. 
This is particularly true in the case that the input signal is the dominant noise in the system and 


thus sets the fluctuation level. In this case the improvement in sensitivity approaches V2. 


3. Radiometer With Dissimilar Amplifiers 


In section 2 the assumption was made that the two amplifiers were identical, except that they 
were allowed to have different effective temperatures. Here the more general problem, in which 
the two amplifiers may differ with respect to other parameters, will be considered. 
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The considerations are limited at first to a case in which the two amplifiers have the same shape 
of gain function but with differing gain amplitudes. If H;(w) is the complex voltage gain of the 
first amplifier, then the complex voltage gain of the second will be taken to have the form 


HAw, t)= [1+ oft) ] Ai). (38) 


where a(t) is a stochastic function. 
If x(t) is the signal resulting from amplifying the signal from the source X (this signal is ((t)) 
with the first amplifier and x(t) is the result of amplifying ((t) with the second amplifier, then a 


straightforward Fourier analysis shows that 
xo t)=[1 + aft) |xi(t), (39) 


and a similar relationship occurs for the signals y;(¢) and y2(t) that arise from the signal m(t) from 
the source Y. This notation is summarized in figure 2. 


The most general form of the autocorrelation function of the sum of two random functions is 





Wa+p = [a(t) + W(t)| [a(t +7) + Wt+7)]. 





= Wa + by + alt) b(t +7) + W(t)a(t +7). 
From this it follows that 


Wear s+ys) (x2 Y2) =Warr2 T Waiyet Wyre ” Wy sus —2x 1 x2Vi1 V2 





—xi(t)yilt + t)ylt)xdlt +7) + yilt)ai(t + rxAt)yelt+7). 


From (39) it follows that 





Ware = X(t) xo(t)xi(t + T)xA(t +7), 





=a t)xX{t+7)[1+a(t)|[1+at+7]. 


= Ware 1+ We + 2a). 


By means of a similar line of reasoning, every term in (41) can be put into a form similar to 


that in (42). Thus, 


Uh 2y+y2r2—Y2) — UK rit+yWri—yyh Tee Wa + 2a), 


l+m+n, 


Vy =X ty, +2, Smoothing 


~ 1 
> Filter 





FIGURE 2. Radiometer circuit with terminology for sec- 
tion 3. 














thes(22—ys) = (We,r, + We,y,)(1 + Ya + 20), 


and in the same way 


Woh are Y2) = (Warr, i Wery )( 1 + a + 2a), 


We,2: = wz 2(1 + Wa + 2a), (46) 


where z!} is the signal that would result if the noise signal n.(t) (the equivalent noise generated in 
the input of the second amplifier) were amplified by the first amplifier. 
When these terms are all collected together, the result is 


Weve =(1+ Wa + 2a){(x?— yi)}t+ We, oF Wi, re Weir, a Wey, a Were, 7 Wey + Were } 
Notice that in this case the second and third terms of (41) are cancelled by the sixth term. 
In order to simplify the analysis, we will consider two cases: (1) a a constant, independent of 
time, and (2) a(t) a random variable with a=0. 


Case I. When ais a constant independent of time, it is elementary to show that 


Woe, =U + a)? { (x3 — y9)? + UWE +) Wee they TY, +o, +, ° (49) 


+1 
217) 


Therefore, both the d-c and the fluctuation terms are multiplied by (1+ a@)?. Since this affects 
both the fluctuations and the deflection by the same amount, the sensitivity is not changed. 

Case II. Assume that a(t) is a random function of time, such that a=0. At balance the term 
(x? — y?)=0, so using a prime to denote quantities evaluated at the balance point, the result is 


Whvg= (1 + ta) {2(0W2, +2 + bee they FH, th (50) 


»1 
-% m4) 


From this the spectral density is 


Qr,v, =(1 + Qa)*{2Qr,*Qz, + Qy,*Qy,) + Qz,*(Qr, +Qy,)+Q,, *(Qz, + Qu)+Q,,*Q:,}- (51) 


Notice that (51) is similar to the result that is obtained in the case of identical amplifiers; the only 
difference is that the present result is the convolution of the original result with (1+ Qa). 
The spectral density of the multiplier output is 


Qu = a*Qy »,, 


16] 





where Q,.., is given by (51). The spectral density at the output of the low pass filter is 
W(f)=G/)QAf), 
= @G(f)(1+ Qa)*{2Qzr,*Qy, + Qy,*Qy,)+Qz,*(Qz, + Qy,) + Qi, *(Qr, + Qy,)+Q.,*Qz,} 
(52) 


Following the same line of reasoning that was used from (11) to (14), the root mean square fluctua- 
tion on the output can be seen to be 


ene = | ; Wfrdf| 


1 
J 


» 


| a2G(f)(1 + Qa)*{2Q0r,*Ox,+ Qu, *Qy,) + Qe, (Oz, + Qy;) 


» 


1 
+Q1 *(Qz, + Qy,)+Q..*Qz,} df . (53) 


To see the effect of the gain fluctuations, consider that the amplifiers have a square bandpass 
of bandwidth B. Then, 


Qr, *Qr, = 207 R2B, 
where R,(w) = H,(w) by (w), 


is the power gain of the amplifier, and 


QaxQr,*Qz, =|" 20?R2BQ o(f—f' df’ =2Q3R2Ba?. (54) 


Thus the term (1+ Q,)*(Q,,*Qz,) can be written as 2Q7R7B(1+ a2). Each term in (53) will produce 
the same effect, so the result is 


Seles v2) 1/2 a9) “ 
W rms = (1 + a?) "2 Wemd 9), (55) 


where w,m40) is the fluctuation level that occurs when there are no gain fluctuations, i.e., when 
a=0. 

Since gain fluctuations increase the fluctuation level of the output but have no effect on the 
average deflection, it is evident that this effect results in a reduction in radiometer sensitivity by 
the factor (1 + a?)!/?. 

While the above expression for the reduction in radiometer sensitivity has been derived only 
in the case of white noise and amplifiers with square bandpass, examination of the integral that 
leads to (54) shows that the result will be qualitatively similar even in a less idealized case. 

The next consideration will be to calculate the sensitivity of a radiometer with differing ampli- 
fier gain functions: (transfer functions). In general these amplifiers could have differing center 
frequencies, and also the shape of the gain-frequency functions could differ. 

The average square output can be expressed as 


w(t)? = | ; G(f)Qudf 
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where G(f) is the power transfer function of the low pass filter. The input to the filter is given by 
U= aVyV2= a(X1 + ¥1 + 21) (X2 — Yo + 22). (2) 
By means of the theorem expressed in (4), the autocorrelation function of u can be put into the form 


a~ "Wy =(xix2 = yiy2)? ait Wr + Wy, * We ) Wr, + Wy, + We.) 





+ X1(t)xX2(t+ 7) X2(t)xi(t+ 7) + yi(t)ye(t+ 7) yo(t)yi(t+7) 





— X2(t)xi(t +7) yi(t)yo(t + 7) — x4(t)xo(t +7) yo(t)yi(t+7). (57) 

Recalling that 

Or =RiQi= Ai AF Q, 
and similar relationships exist for the other variables, the spectral density of u becomes 
a? Qu =(X1X2— V1 ¥2)* Hf) +(RiQr+ RiQmt Ri Qn rl R2Qi+ R2Qm+ RQ) 
+ Re {{(H, HF Q))*(H2H¥ Q)| + (Hi HS Qn)*(H2H* Qn)| 
— [(H2H% Q))*(H1 HF Qn)| — (41 H3 Q)*(H2H* Qn) |}. (58) 
The form of the last term arises from the theory of cross-spectra [Goodman, 1957; Korn and Korn, 
1961}. 

The complete expression for the spectral density of the multiplier output is obtained by 
evaluating (58). From the resulting expression the fluctuation amplitude of the radiometer can 
be obtained by following the same procedure used previously. 

Next it is necessary to calculate the average deflection resulting from a particular combination 
of input signals. This can be obtained by evaluating (2) of the simple theory in this more general 


case. 


First, consider 





U = AV V2 = A(X + ¥1 +21 )(X2— V2 +22). (2) 


“ 


Of the six random variables that appear in the right-hand term of (2), the only correlations 
that are not zero are the partial correlations of x; with x2 and of y; with yz. Therefore this expres- 
sion reduces to 


u = a(x1X2 - Y1¥2). (59) 
Equation (59) must be evaluated in terms of the transfer functions of the two amplifiers, 


H,(f)+H2(f), and the spectral densities of the noise outputs of the two noise sources. The Fourier 
transform of the random function / is 


Sif) =| l(t)e-—27¥'dt. (60) 


In the theory of cross-correlations, it is shown that [Korn and Korn, 1961] the time average of 
the cross correlation of x; and x2 can be expressed as 


T 


(xix) = | Re| Lim 7 S288) |d/= Lim 1. | (Sz,S4,+ S3Sx)df=3 | (H\H¥+H¥H)Qdf, (61) 
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where the symbol ( ) is used to indicate a time average. The final form of (61) results from the 
fact that S; =H,S), etc. Also, due to the fact that the random processes are ergodic, the time 
average in (61) may be replaced by an ensemble average. 

In the same way 


Ty | (MHS + H2HF)Qmndf. 


Therefore, 
u= a[x1xX2— V1 y2]= - | (H,H* + H2H*)(Qi—Qn)df. (63) 


The above expression leads to the following conclusions: (1) In case the two sources have the 
same form of spectral density over the amplifier bandwidth, balance occurs independent of the 
amplifier response curves. Thus, for white noise, balance occurs when 7;=T7,,. (2) In case Q; 
and Q,, have different forms, balance occurs where 


[ ens + Hatt Qrdf= | HLHE+ H2HF)Qndf. (64) 


An expression for Au, the change in u due to a deviation from balance, is required. In the 
case where both signals consist of white noise, 


| —" 
Q1— Qm=4MTi— Tn). (65) 
Next, the assumption is made that 
Tn=1,— AT, 
I" a 
then. Qi— Om=Z ATi —(T— AT)|= j RAT. (66) 
In this case 


Au=5 KAT | (HHS + H2H*)df. 


If the smoothing filter has a response function at zero frequency given by Go, then 


Aw=G12Au a GhkAT | (H,H3 + HH df. (68) 


Again the assumption is made that the minimum detectable change in deflection is a change 
equal to the rms value of the fluctuations. The minimum detectable value of AT is evaluated by 
equating the right-hand side of (68) and the square root of w?. Using the fact that 


=the ‘fe 
Wrms—(w*) = | wi fds| ; 


the result is 


x 1/2 
| CHOUNas| =" GY2kAT | (H,H¥ + HeH* df. 
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When this solved for AT, the result is 


% 1/2 
| [ enewnar 
AT= a 
cyt | (H,H* + H2H*)df 





a 


8 


Equation (70) is the general expression for the sensitivity of a correlation radiometer with 
amplifiers with differing complex gain functions. In order to interpret this result and obtain a 
feeling for the importance of the various parameters, it is convenient to look at two examples. 
First the case of two amplifiers with square gain functions centered on the same frequency fo; the 
first amplifier with bandwidth B, and gain R? and the second with bandwidth By, and gain R$ will be 
considered. To be specific, it will be assumed that B; > By. It will be assumed that no phase 
shifts occur in these amplifiers. 

In order to evaluate (70) it is necessary to evaluate the various terms of (58), recalling that the 
numerator of (70) contains only fluctuation terms; i.e., any d-c contributions must be ignored. It 
is reasonable to evaluate the convolution integrals that occur with the restriction that fis small 
enough that B, is enclosed within B;. When these evaluations are carried out with the assumptions 
that the radiometer is balanced and also that the low pass filter has a square passband of width b, 
the result is 


[numerator of (70)|= aGl/?b"*(B.R°R8)'?k(T) + T,). 


When this is compared to the similar result from the simple theory, for example (14), it is evident 
that the bandwidth that sets the fluctuation level is B., the narrower of the two amplifier bandwidths. 
Also, the amplifier gains appear as the product R°R9. 

The denominator of (70) must also be evaluated under the same conditions. This involves 
evaluating the integral 


[ unne+ nampa (71) 


where H, and H» are the discontinuous real functions as defined above. When this is carried 
out, the result is 


as pare ' 
[denominator of (70)|=5 Gi? kK(4R¢RSB3). 


These results lead to 


1/2 
aT=2(2 (T,+T,). (72) 


Thus, in this case the same sort of expression is obtained as in the case of the simple theory 
except that the bandwidth of the narrowest amplifier sets the sensitivity. 

Finally the case of two amplifiers with square bandpass with identical bandwidth but tuned to 
slightly different center frequencies will be considered. Again phase shifts will be neglected, the 
complex gains will be treated as real positive functions, and the radiometer will be assumed to be 
balanced. The resulting sensitivity is 


ij 1/2 
AT=2 fj (Ti+-T 


where By is the “overlap bandwidth” or the frequency interval common to both amplifiers. 
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4. Effect of Variation in Time Delay and Phase in the 
Amplifiers of a Correlation Type Radiometer 


Here it will be assumed that a differential variation in time delay can occur to the signals in 
the two amplifiers. Thus, at any instant the multiplier is comparing yoltages corresponding to 
two different instants of time. Therefore, the multiplier output can be expressed as 


u(t)=v,(t)ve(t+7), (74) 


where 7 is the time delay difference in the two signal channels. In order to simplify this analysis, 
the assumption will be made that the two radiometer channels are identical except for the dif- 
ference in time delay. Thus the voltages that are multiplied together are 


V(t) = x(t) + Vib) + 2Z(F), 


DAL TT) = ME TT) HET TFT) TT 2akt TT). (75) 


Equation (7) shows that the fluctuation level of the radiometer depends only on convolutions 
of power spectra. Since these are not a function of 7, the fluctuation level does not vary as a result 
of a differential time delay. Therefore, the change in sensitivity will be related to the change in 
average deflection caused by the time delay. 


The change in average deflection can be obtained from 


D=([d)+v)+210] + r)— y(tt+7)+ z2(t+7)] 


(talt+7) — MOMt+ 7) =WA7) — W,/7). (76) 


The Wiener-Khintchine theorem allows this to be put into the form 


u(t, T) = (0,—O eas (77) 


In order to be able to evaluate the above integral. it is convenient to consider that both sources 
generate white noise and the amplifiers have square bandpass of width B centered on fo. Thus. 


in this case, 


i 3 + B/2 


{ j H/2 ‘ Jin 
m7) =(Q,—Q,) | eld f + | e" rdf| / 
LJ —fo—B/2 fo-B/2 


si a77B 
= 2B(\Q0,,—Q,) cos 27Tfo ae 


This can be put in more convenient form by noting that 
u(0) = 2B(Q, — Q,). 


sin wrB 
77B 


u(t) = uO) cos 277fo 
For the usual case in which B < fy and 7 is on the order of 1/fo, this can be approximated by 


u(t) = ul) cos 277fo. 
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Thus, in the case of white noise sources and square bandpass amplifiers, time delay variations 
will not influence the balance point. They will affect sensitivity, however. For maximum sen- 
sitivity, it is important to keep the difference in time delay in the two amplifiers small enough 
that cos 277fo is approximately unity. 

A somewhat related question is the effect of a differential phase change in the two amplifiers. 
Again assuming white noise sources and amplifiers that are identical except for the phase shift, 
the average deflection can be computed by means of (63). With the assumption that 


Rs = R,e??, 


this equation becomes 


a 
u(d) —s (Q:—Qm)R (2 cos @). 
This can again be expressed in the form 


u(d@) = u(0) cos @. 


5. Some Effects of Using an Imperfect Multiplier 


The simplest form of multiplier to multiply two microwave signals together is probably that 
shown in figure 3. The input signals (v; and v2) are applied to two opposite arms of a hybrid junc- 
tion. The other pair of opposite arms are terminated by square law detectors. The outputs from 
these detectors are applied to a difference amplifier. In case the hybrid Tis matched, the detectors 
have the same sensitivity and their reflection coefficient is zero: and, if they are assumed to have 
a perfect square law envelope response, the output voltage z will be 


(" + ) (* — sy" 
2=],—- =a] eh "I. 
te V2 Vz 


where a is the coefficient of proportionality in the multiplier law. Thus, 


z=2avI>. (84) 


In an actual multiplier constructed on this principle, the ideal conditions suggested above will 
be only approximated. Thus, the output will not be an accurate multiplication of the input signals. 
A scattering matrix analysis of the multiplier junction indicates that if the junction is matched and 
perfect, the wave arriving at the detectors is equal to that assumed in the ideal case as long as 
either pair of leads are terminated in matched loads. Since the detectors (if bolometers) can be 
quite accurately matched, this will probably not be the phenomenon that sets the multiplier 
accuracy. This accuracy can depend on the accuracy of the “‘square law” of the detectors, 
however. 





a FIGURE 3. Microwave multiplier. 
+——— Z 


* 





Differential 
Amplifier 
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In order to investigate the effect of the detector law on the multiplier accuracy, assume that 
hoth detectors follow the same law. If one detector is more sensitive than the other, the incoming 
signal can be attenuated, so that, referred to the input of the attenuator, the two detectors are 
assumed to have the same law with the same numerical coefficients. Thus, the voltages from the 
detectors can be expressed as power series 


(85b) 


Then the multiplier output is 


These terms may be expanded with the binomial theorem: the result is 


ds ° , a4 . 2 
Oe pe a EE fp a SS See ye beet 3 4 1.9) + ° - 
(2vv2) + 55 (Buje + vy) + S75 (Aviv. + 4uiv)) + . . . - (87) 


» 


92/2 


The smoothing filter provides an estimate of the average of the multiplier output: thus it is neces- 
sary to compute z._ This is, 


z= 2[ ao( v2) T a Viv2 T v\v3) “ee 








WaAx + y + 21)(x — y+ 22) t+ ag{(x + y + 21) a — y+ 22) + (x + y+ 21)\(x— yt 22)5} 


Balx* — y*) + 2agx*— y*) + 3adx? — y*) (22 +23)+. . .]. (88) 


The result is that if x‘—j4=0 at the same conditions that x*7— y*=0, the balance condition 
is not affected by any term up to the fourth in the crystal law expansion. This would be the case 
if both x and y are signals with the same statistical properties: for example, if both are thermal 
noise sources, they both possess Gaussian statistics. 

However, in the case that the two signals have different statistical properties, (x*— y*) would 
not necessarily equal zero when (x2— y°) does. Thus, in this case, it is important to select de- 
tectors that are accurately square law. 


6. Conclusions 


The sensitivity of a correlation radiometer has been computed under a variety of conditions. 
In section 2 a very simplified calculation is used to derive the usual expression for the sensitivity 
of the radiometer. This is followed by a calculation of the sensitivity of a radiometer of the type 
suggested by Allred (1962) in which the unknown noise signal is balanced against a sinusoidal 
reference signal. By performing these two calculations in an analogous manner, it is particularly 
easy to compare the sensitivities that result. In the case that the output fluctuations are pre- 
dominately due to the input noise signals, Allred’s radiometer is more sensitive than the conven- 
tional circuit that compares two noise signals. 
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The first part of section 3 demonstrates that it is not necessary for the two amplifiers to have 
the same gain. The calculation shows that the product of the gains is the parameter that deter- 
mines the output amplitude, instead of the individual amplifier gains. Next, it is shown that gain 
fluctuations produce a decrease in sensitivity. The last part of this section is concerned with the 
effect of using two dissimilar amplifiers. A general expression is obtained for the sensitivity of a 
radiometer with amplifiers with arbitrary gain functions. This expression is evaluated for the case 
of amplifiers with square bandpass and no phase shifts. It is shown that under these conditions, 
whether the two amplifiers have the same bandpass and different center frequencies or whether 
they have the same center frequency and different bandwidths, the sensitivity is determined by 
the “overlap” bandwidth. 

Section 4 deals with the effect of a differential time delay or differential phase shift in the two 
radiometer channels. It is shown that both of these effects result in a decrease of sensitivity. 

Finally in section 5 it is shown that, if the multiplier carries out the multiplication operation 
by forming the difference of the squares of the sum and difference of the two input signals, errors 
can result if the “square law” elements do not have a perfect square law response and if the two 
input voltages have differing statistical distributions. Thus, the radiometer proposed by Allred, 
in which the comparison signal is sinusoidal, places a more stringent requirement on the multi- 
plier than does the more usual correlation radiometer, in which both input signals are Gaussian 
noise. 


The author gratefully acknowledges many helpful discussions with D. F. Wait and G. F. Engen. 
M. M. Siddiqui’s comments on the manuscript proved to be very helpful. 


7. References 


Allred, C. M. (1962), A precision noise spectral density comparator, J. Res. NBS 66C (Engr. and Instr.), No. 4, 323-330. 

Blum, E. J. (1959), Sensibilité des radioteléscopes et recepteurs a correlation, Annales d’ Astrophysique 22, No. 2, 140-163. 

Colvin, R. S. (1959), Faint signal limitations of radiometers, Inst. of Radio Engr., WESCON Convention Record 3, Pt. 8, 
52-58. 

Colvin, R. S. (1961), A study of radio-astronomy receivers; Stanford Radio Astronomy Institute Publication No. 18A (Stanford 
University, Stanford, Calif.). 

Goodman, N. R. (Mar. 1957), On the joint estimation of the spectra, cospectrum and quadrature spectrum of a two-dimen- 
sional stationary Gaussian process, Engineering Statistics Laboratory Scientific paper No. 10, New York University, p. 23. 

Korn, G. A. and T. M. Korn (1961), Mathematical Handbook for Scientists and Engineers (McGraw-Hill Book Co., Inc., 
New York, N.Y.). 


8. Appendix — Evaluation of Convolution Integrals 


Assume a random signal with white spectral density amplified by an amplifier with square 
passband as shown in figure A-1. 

If the spectral density of the input signal is Qo as measured, the amplitude of Q; is RQo/2, 
where R is the power gain of the amplifier. The division by 2 occurs because the spectral density 
is assumed to be split equally between the positive and negative frequency regions. Then, 


0,*Q,= [ Of FOP df’. 


J 
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FIGURE A-1. Square spectral density. FIGURE A-2. Functions appearing in integrand of con- 
volution integral: 
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The value of this integral for small values of the parameter fis required. The term Q;(f—/’) has 
the same shape as Q,(/') except that it is displaced an amount f and inverted on the frequency 
axis, as shown in figure A—-2. Since the integrand is the product of these two functions, the value 
of the integral is the product of the amplitudes of these functions times the frequency interval 
over which the product is nonzero. Thus, 


Q2*Q2=(RS) 2B-/) 


for f < B. 
In this paper, the value of this integral is required under the condition that f<B. To this 
approximation, the result is 


2 


Q:+0.=28(R&) (A-2) 


In general, the convolutions of the other spectral densities appearing in this paper are evalu- 
ated in the same way. The reason that the various convolutions of power spectra that are used 
have differing numerical coefficients is that the hybrid junction at the input of the radiometer 
divides the input powers; however, the noise powers appearing due to the effective input tempera- 
tures of the amplifiers do not undergo this power division. 
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Research on crystal growth and characterization at the 
National Bureau of Standards June 1966, ed. H. C. Allen, Jr., 
Tech. 293 (Aug. 26, 1966), 25 cents. 


The National Bureau of Standards with partial support from the 
Advanced Research Projects Agency of the Department of Defense 
is continuing a wide program of studies involving crystalline mate- 
rials. These include investigation of methods and theory of growth, 
study of detection and effects of defects, determination of physical 
properties, refinement of chemical analysis, and determination of 
stability relations and atomic structure. The types of materials 
range from organic compounds, through metals, and inorganic salts 
to refractory oxides. This report summarizes 
projects wholly or partially supported by ARPA. 


progress in those 


Key Words: Crystal growth, crystal characterization, and crystalline 
materials. 


A sensitive recording NMR ultrasonic spectrometer, L. 
James, Tech. Note 344 (Sept. 7, 1966), 20 cents. 

Instrumentation is described for automatically obtaining a con- 
tinuous recording of ultrasonic absorption lines. 


W. 


Ultrasonic power 
is added to the crystal lattice by a transducer mounted on a single 
crystal of the sample. The interaction of the ultrasonics with the 
nuclear spin system is recorded by using a pulse NMR system in 
conjunction with a boxcar integrator. This technique is more 
sensitive and provides more easily interpretable results than earlier 
systems. 


Standard Reference Materials: Viscosity of a standard lead- 
silica glass, A. Napolitano and E. G. Hawkins, Mise. Publ. 260-11 
(Nov. 7, 1966), 25 cents. 

The viscosity of a lead-silica glass has been measured at the National 
Bureau of Standards and seven other laboratories. Determinations 
were made in the range of 10? to 10! poises (1350-400 °C). Meas- 
urements were made by the rotating cylinder, restrained sphere, 
fiber-elongation, and beam-bending methods. The results have 
been critically evaluated and the glass has been issued as Standard 


Reference Material No. 711. 


Key Words: Beam bending, fiber elongation. glass, glass standard. 
glass viscosity, lead-silica glass, restrained sphere, rotating cylinders, 
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Thermal conductivity of selected materials, R. Powell, 
C. Y. Ho, and P. E. Liley, NSRDS-NBS 8 (Nov. 25, 1966), $1.00. 
The work presented in this report consists of the critical evaluation 
and analysis of the available thermal conductivity data on eleven 
metals and nine nonmetals for the solid state, on seven fluids for 
both the liquid and gaseous states and on two for the liquid state only. 
The materials studied were selected primarily for their potential 
applicability as reference standards or because of their technical 
importance. The temperature range for which values are given 
often exceeds that for which these values are known with a high 
degree of certainty. 


Key Words: Thermal conductivity, solids, metals, liquids. gases. 
Notes on the state-of-the-art of benefit-cost analysis as related 


to transportation systems, J. D. Crumlish, Tech. Note 294 (Nov. 1, 
1966), 30 cents. 
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This review of benefit-cost analysis as a tool for evaluating alternative 
courses of action describes the technique, discusses a number of 
benefit-cost studies, and indicates the difficulties inherent in this 
area of applied economics. The author concentrates on the applica- 
tion of the technique to large scale transport problems, reviews the 
literature and indicates in his conclusions where the technique can 
be helpful and where there is little chance for its success. 

An accompanying matrix of benefit-cost studies and a commentary 
thereon is supplied by Marsha Geier, an NBS economist. Miss 
Geier’s literature search failed to produce any analytic methods 
which were comprehensive, theoretically justifiable, operational or 
significant. This finding tended to support the views of the author. 
(Introduction by Alan J. Goldman, Applied Mathematics Division, 
Institute for Basic Standards, National Bureau of Standards.) 


Key Words:’ Benefit-cost analysis: transportation economics; trans- 
port systems; systems analysis: state of the art. 


Calculation of the heating value of a sample of high purity 
methane for use as a reference material, G. T. Armstrong, 
Tech. Note 299 (Dec. 15, 1966), 25 cents. 

The heat of combustion of CH, has been recalculated in kJ (mol)~', 
Btu (mol)~', Btu (cu ft)~! (dry basis) and Btu (std. cu ft)~ ! (saturated 
basis), using the best available experimental determinations of the 
heat of combustion and other measured quantities and the most 
recent generally accepted physical constants and defined physical 
units. The calculations are outlined in detail. The resulting quan- 
tities are applied to calculation of the heat of combustion of a 
reference sample of CH, submitted for analysis of composition and 
certification of heating value by the Institute of Gas Technology. 


A standard for accurate phase-angle measurements at audio 
frequencies, W. W. Scott, Jr.. Tech. Note 347 (Oct. 14, 1966), 
25 cents. 

\ method is described for the measurement of phase angles from 
0° to 360° (except for a band + 10° about the 180° point) with an 
uncertainty at audio frequencies of +0.05° to +0.01°, depending on 
the phase-angle. The method utilizes readily available equipment 
which may be assembled and connected to serve as astandard. The 
principal advantage of this standard is in its potentially broad fre- 
quency-response; namely, 10 Hz to 10 kHz or higher. It makes use 
of inductive voltage dividers and thermal voltage converters to com- 
pare voltage magnitudes, which are related to phase-angle by the 
law of cosines. The thermal voltage converters are compared with 
each other before each phase-angle measurement. Their outputs 
are connected so that short-time power supply variations do not 
affect the measurements. The repeatability of measurements with 
this standard at a phase-angle of 60° is within a range of 0.001°. 
A model has been built which shows good measurement precision, 
but it has not been thoroughly tested or studied at many frequencies 
and phase-angles for sources of error. This report serves as a 
record of the development work completed on this standard to date. 


Key Words: Audio frequency. inductive voltage divider, measure- 
ment. phase-angle, thermal voltage converter, standard. 


A new near-zone electric field-strength meter, F. M. Greene, 
(Proc. 8th IEEE Symp. Electromagnetic Compatibility, San Fran- 
cisco, Calif., July 11-13, 1966), Symposium Digest New Ideas for 
Electromagnetic Compatibility, 1966. 

The National Bureau of Standards has recently completed the 
development of prototype instrumentation for measuring the electric- 
field components of complex, high-level, near-zone electromagnetic 





fields from 0.1 to 1000 volts per meter, at frequencies from 150 kHz 
to 30 MHz with a present uncertainty of less than +2 dB. The 
successful design of the NBS meters is based on the use of a novel 
form of telemetry, employing a completely non-metallic electrical 
transmission line, which apparently has not been fully exploited 
heretofore. This avoids the perturbing effects on the field being 
measured, usually caused by field-strength meters employing metal- 
lic ore transmission lines. 


Key Words: Device, electro-explosive ordnance —field, near-zone 
electromagnetic —hazards, electromagnetic radiation—line, non- 
metallic electrical transmission—meter, electric field-strength— 
telemetry, novel form. 


New standard for the calibration of microphones, W. Koidan, 
Mag. Std. 37, No. 5, 141-144 (May 1966). 

This article discusses the background, contents and use of a new 
American Standard, Calibration of Microphones, $1.10—1966, pub- 
lished by the American Standards Association, Incorporated. 


Key Words: Calibration, microphones, Standard. 


Resistance to weathering, H. T. Arni, ASTM Significance of Tests 
of Concrete and Concrete Making Materials, Am. Soc. Testing Mater. 
Spec. Tech. Publ. 169A, pp. 261-274 (1966). 

Deterioration of concrete due to physical processes produced by 
exposure to natural weathering is discussed. The most recent 
theories on the mechanisms of freezing-and-thawing damage are 
presented, and freezing-and-thawing tests now specified and _ pro- 
posed are discussed. The effects of de-icing chemicals and of 
weathering processes other than freezing and thawing are also 
considered. 


Key Words: Concrete, de-icing scaling, durability, freezing-and 
thawing, weathering. 


Survey of infrared measurement techniques and computa- 
tional methods in radiant heat transfer, S. T. Dunn, J. C. Rich- 
mond, and J. F. Parmer, J. Spacecraft Rockets 3, No. 7, 961-975 
(July 1966). 

This paper is an analysis of 133 replies to a questionnaire sent to all 
known workers in the field of thermal radiation property measure- 
ments and radiant heat transfer computations. The questions 
were designed to obtain information on the current state of the art 
in these fields. In addition, the authors have supplemented the 
information obtained from the questionnaire to give a broad overall 
picture of the current status and future needs in each area covered. 
4 brief summary of Russian literature in the field of thermal radia- 
tion property measurements is also included. 


Key Words: Computational methods, emissivity, emittance measure- 
ments, radiant heat transfer, reflectance measurements, thermal 
radiation. 


Optical methods of temperature measurement, G. A. Horn- 
beck, J. Appl. Opt. Lead Article 5, No. 2, 179 (Feb. 1966). 

The history of the laws of thermal radiation is briefly reviewed, and 
the theoretical relations basic to some of the optical methods of 
temperature measurement are given. Comments on a number of 
optical methods are presented, as well as a selected reference with 
each method where more detail is available. In conclusion, several 
problems, in areas of uncertainty in the field of optical methods of 
temperature measurement, are suggested in the belief that a solution 
to these problems will greatly enhance the technical value of this 


field. 


Key Words: Boltzmann, Draper, Jeans, Kirchotf, optical, Planck, 
pyrometry, radiance, Rayleigh, Stefan, temperature, Wien. 


Infrasonic waves from aurorae, K. Maeda and J. M. Young, 
Nature 207. No. 4994, 279-281 (July 1965). 

The effects of atmospheric thermal structure on the free-field propa- 
gation of auroral infrasound are discussed. Experimental data 
illustrating dispersion and near-sinusoidal oscillations near the 
Bruunt frequency are shown. The effects of the large positive 
lapse rate layers, especially the one in the thermosphere, on atmos- 
pheric resonances are given particular attention. 


Magnetic core permeability measurement techniques, R. D. 
Harrington and A. L. Rasmussen, Proc. Magnetic Core Conf. 7, 
11-24 (1965). 

Techniques for measuring the initial complex permeability of mag- 
netic cores in the frequency range from 1 kHz to 50 MHz are briefly 
discussed. Emphasis is directed toward studies recently carried 
out at the National Bureau of Standards. 

At frequencies below 100 kHz, either a demountable coil or a spe- 
cially designed Maxwell bridge is used for obtaining core data. The 
frequency range from 100 kHz to about 50 MHz is most conveniently 
covered through the use of the radio frequency permeameter. A 
brief survey of the development of the permeameter with a descrip- 
tion of some recent measurements indicating that this method may 
be used to frequencies even greater than 50 MHz is given. 

It has been found that toroidal coils are preferable to the perme- 
ameter for obtaining accurate measurements of the temperature 
coefficient of permeability. A brief description is given of a method 
for evaluating the inductance and resistance changes of such coils 
using a bridge as a null indicating device. The technique is capable 
of measuring temperature coefficients of permeability of the order 
of one 10-® per degree centigrade. 


An examination of the effects of heat transfer and compres- 
sible flow on the performance of laminar flowmeters, F. W. 
Ruegg and H. H. Allion, Proc. Flow Measurement Conf., Am. Soc. 
Mech. Eng. Fluid Meters Golden Anniversary, Pittsburgh, Pa., 
Sept. 26-28, 1966, pp. 253-273 (1966). 

Laminar meters are now widely used for measurement of gas flow, 
and high accuracy can be achieved with proper calibration and use. 
It has been found necessary to modify Poiseuille’s law for laminar 
flow to explain the performance of meters over the wide range of 
flow conditions encountered. A one-dimensional flow analysis 
is used to derive the effects of meter shape, heat transfer, and com- 
pressibility of the gas on meter performance. A relationship be- 
tween the meter flow coefficient and the similarity parameters of 
Reynolds number, Prandtl number, Mach number, and Knudsen 
number is given in a form that is convenient for both the calibrator 
and user of the meter. The small but significant effects predicted 
are compared with experimental results. 


Key Words: Compressible flow, flowmeter, gas meter, heat transfer, 
laminar flow, meter calibration. 

Roofing research, T. H. 
(Sept.-Oct. 1966). 

The Building Research Division of the National Bureau of Standards 
has for many years studied roof systems, constituent materials and 
components. Highlights of the roofing research in the last five 
years and a brief résumé of the results of this research are given. 
Description of some of the newly developed roofing systems and the 
present status of roofing research are also discussed. 


Boone, Military Engr. 385, 344-345 


Key Words: Liquid-applied 
prefabricated sheets: 
heating. 


membrane movements; 
radiative cooling: roofing research: solar 


coatings: 


A new guide for photographers, C. S. McCamy, Mag. Stds. 37, 
No. 8, 223-224 (Aug. 1966). 

The new “American Standard Photographic Exposure Guide” 
presents and relates data which provide means of calculating camera 
settings for daylight and moonlight exposure of black-and-white 
and color films. The Appendix gives guidance for photographing 
sunsets, clouds, the moon, artificial satellites, rainbows, cities at 
night, fireworks, lightning, aurora, stars, and eclipses. The standard 
is based on extensive picture tests. 


Key Words: Camera, crystal, diffraction, dislocation, topography, 
X-ray. 


A simple non-scanning camera for x-ray diffraction contrast 
topography, H. P. Layer and R. D. Deslattes, J. Appl, Phys. 37, 
No. 9, 3631-3632 (Aug. 1966). 

A new x-ray topographical camera has been developed that is appli- 
cable to both Laue and Bragg geometries for normal and anomalous 
transmission, that results in a significant reduction of experimental 
complexity along with quite acceptable exposure times. 


Key Words: Camera, diffraction, dislocation, topography, x-ray. 
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Studies of some exploding wire light sources, E. C. Cassidy 
and S. Abramowitz, J. SMPTE 75, No. 8, 735-737 (Aug. 1966). 
Continuous and time-resolved measurements of the spectral dis- 
tribution of light emitted by various exploding wire systems have 
been obtained by use of a high-speed drum camera and a rotating 
shutter, respectively. Results from experiments with several sys- 
tems show the effects of environment, pressure, energy, and wire 
material on the spectrum. Various intermediate species, produced 
by the explosion, are determined spectroscopically. 


Key Words: Al0, AlH, aluminum, atomic spectra, electrical dis- 
charges, exploding wires, light sources, time-resolved spectroscopy, 
titanium TiO. 

A device for use in applying strain gages to cylindrical speci- 
mens, T. W. Butler and R. L. Bloss, Experimental Mechanics 6, 
No. 10, 528 (Oct. 1966). 

A simple device is described which applies a uniform pressure to 
the exterior of a cylindrical specimen. This device has been used 
in connection with a regulated air supply to apply pressure while 
curing strain gage installations. 


Key Words: 


strain gage. 


Adhesive, application, curing, cylinder, installation, 


Allion, 
, 1965), 


Calibration and use of laminar flowmeters, H. ped 
(Proc. 25th Annual Gas Measurement Short Course, Aug. 2 
West Virginia Univ. Tech. Bull. 77, pp. 34-43 (1965). 
The laminar meters discussed are made of one or 
passages each having a depth or diameter small in relation to its 
length. Resistance to flow is largely viscous and may be analyzed 
in a fairly simple fashion in the incompressible case. For gas flow 
in these meters, treatment of the effects of deviations from incom- 
pressible theory is quite complex; however, errors resulting from 
ignoring these deviations when predicting meter response at other 
than calibration conditions are apparently less than the measure- 
ment errors over a fairly wide range of conditions. In some appli- 
cations linear response and independence from pressure level may 
be approximated closely. 


more uniform 


Key Words: Accuracy, 
meter and effects, 
Mach number, 


calibration, compressible flow, differential 
fluid meter, incompressible flow, laminar flow, 
Reynolds number, viscous flow, volumetric flow. 


Embrittlement of high strength AISI 4340 steel in boiling 


NaCl solution, H. Logan and J. 
No. 9, 265-269 (Sept. 1966). 

‘he resistance of high strength AISI 4340 steel tubing, quenched 
in oil and tempered at 425 °F, to corrosive media has been investi- 
gated. The exterior surface of the specimen was exposed to boiling 
NaCl, the interior cavity was evacuated and the specimen was sub- 
jected to a tensile Specimens failed, with little evidence 
that cracking was developing, after several hours exposure and after 
hydrogen had penetrated through the specimen wall. The appear- 
ance of the initial fracture and its mode of development indicate 
that it originated within the specimen wall. The data indicate that 
failure resulted from hydrogen embrittlement. 


M. Wehrung, Corrosion 22, 


stress. 


Key Words: Boiling sodium chloride solution, 
strength steel, hydrogen embrittlement, 
initiation, stress-corrosion. 


delayed failure, high 
hydrogen diffusion, internal 


Some effects of notch geometry on the tensile behavior of 
annealed type 310 stainless steel, ¢ Geil, J. Materials 1, 
No. 3, 583-608 (Sept. 1966). 

A study was made of the influence of notch geometry (angle, depth 
and root radius) on the tensile deformation and fracture charac- 
teristics of circumferentially notched cylindrical tensile specimens 
of annealed 310 stainless steel at ambient temperatures. Speci- 
mens with a root radius of 0.05 inch or less and theoretical stress 
concentration factors of 1.6 or greater, regardless of notch angle 
or depth, exhibited cracking at small plastic strains. These strains 
were generally much smaller than the strains at maximum load. 
The cracking slowly increased in depth in a discontinuous manner 
with increase in applied stress and extension of the specimen until 
a condition of plastic instability was reached, resulting in a final 
rapid fracture of the specimen. Notched specimens with a root 
radius of 0.10 inch or greater, and theoretical stress concentration 


factors less than 1.6 did not exhibit any cracking at the root of the 
notch prior to the final rapid fracture. 


notch 
tensile defor- 


Key Words: Fracture, initial cracking, multiaxial stresses, 
geometry, 310 stainless steel, stress concentrations, 
mation, workhardening. 


A nonmagnetic laboratory for the National Bureau of Stand- 
ards, F. K. Harris, IEEE Spectrum 3, No. 11, 85-87 (Nov. 1966). 
Plans are discussed for a nonmagnetic laboratory to be built on the 
Gaithersburg campus of the National Bureau of Standards for 
absolute electrical measurements and other purposes. The opera- 
tion of the NBS current balances and the proton gyromagnetic ratio 
apparatus is briefly reviewed, and a short history of the national 
electrical units since the creation of the Bureau is included. 


Key Words: Absolute ampere determination, calculable inductor, 
current balance, national electrical units, nonmagnetic facility, 
pellat dynamometer, proton precession frequency. 


A scattering chamber for use with cooled large area lithium- 
compensated silicon radiation detectors, W. R. Dodge, J. A. 
Coleman and S. R. Domen, Nucl. Instr. Methods 42, 181-187 
(July 1966). 

A scattering chamber designed to be used with large area lithium- 
compensated silicon radiation detectors has been constructed... The 
chamber has provision for cooling a circular array of detectors to 
77 °K to achieve short charge collection times and good resolution. 
The apparatus was designed to be used with the NBS 180-MeV 
synchrotron in photoproton experiments in the 2 to 35 MeV proton 
energy interval. 


Key Words: Detector telescope, large area detectors, photonuclear 
physics apparatus, scattering chamber, semiconductor detectors. 


Lithium-compensated silicon focal plane detectors for elec- 
tron scattering spectrometers, W. R. Dodge, J. A. Coleman, S. R. 
Domen and J. K. Whittaker, Rev. Sci. Instr. 37, No. 9, 1151-1159 
(Sept. 1966). 


In order to utilize the inherent resolving power of a double-focusing, 
169.8°, 76.2-cm radius-of-curvature magnetic spectrometer being 
constructed at the National Bureau of Standards for elastic- and 
inelastic-electron nuclear scattering studies, focal plane counters 
were needed whose dimension along the direction of momentum 
dispersion was l-mm. This paper describes the fabrication tech- 
niques and performance of lithium-compensated silicon 1 X 1 X 60- 
mm radiation detectors which seem well suited for use as focal plane 
detectors for the above spectrometer. Differential and integral 
pulse height distributions for 10-40 MeV electrons incident on the 
counters are presented as a function of the temperature and bias 
voltage. Effects due to the finite width of the counters are dis- 
A tunnel diode discriminator which provides the desired 
binary information to the NBS data processing system is described. 


cussed. 


Key Words: Electron scattering, focal plane detectors, lithium- 
compensated silicon, low temperature radiation detectors, photo- 
nuclear physics apparatus, rectangular counters, spectrometer. 


Measurement of the mean energy required to create an 
electron-hole pair in silicon between 6 and 77 °K, W. R. 
Dodge, R. Domen, T. Leedy, and D. M. Skopik, Phys. Rev. 
Letters 17, No. 12, 653-655 (Sept. 19, 1966). 

Relative values of the mean energy required to create an electron- 
hole pair with 7°7Bi conversion electrons between 6 and 77 °K have 
been measured. The data show that € exhibits no first order varia- 
tion with temperature from 6 to 77 °K. Previous measurements had 
yielded values as high as 5.22+ .02 eV/electron-hole pair at 20 °K. 


Key Words: Energy require to create an electron-hole pair, low 
temperature, recombination, Shockley theory of €, silicon radiation 
detectors, trapping. 


Multiple diffraction in the Weissenberg methods, A. Santoro 
and M. Zocchi, Acta Cryst. 21, Pt. 3, 293-297 (Sept. 1966). 

The conditions for multiple diffraction due to symmetry have been 
derived for the Weissenberg geometry, for the various crystal sys- 
tems and for the most commonly used rotation axes. 
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It is shown that the equal-cone method is the most appropriate in 
intensity measurements. 


Key Words: Crystal structure, diffraction, multiple diffraction, 
Weissenberg method. 
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